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Abstract 

Main arguments on the Dark Matter particle direct detection approach are 
addressed on the basis of the work and of the results of the ~ 100 kg highly 
radiopure Nal(Tl) DAMA experiment (DAMA/Nal), which has been operative 
at the Gran Sasso National Laboratory of the I. N.F.N, for more than one decade, 
including the preparation. The effectiveness of the WIMP model independent 
annual modulation signature is pointed out by discussing the results obtained 
over 7 annual cycles (107731 kg • day total exposure); the WIMP presence in the 
galactic halo is strongly supported at 6.3 a C.L. The complexity of the corollary 
model dependent quests for a candidate particle is also addressed and several of 
the many possible scenarios are examined. 
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1 The physical problem 

1.1 Evidence for Dark Matter in the Universe 

The first evidence that much more than the visible matter should fill the Universe 
dates back to 1933 when F. Zwicky measured the dispersion velocity in the Coma 
galaxies [1]. This was soon after confirmed by S. Smith studying the Virgo cluster [2]. 
Nevertheless, only about 50 years later the fact that Dark Matter should be present 
in large amount in our Universe finally reached a wide consensus. 
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Particular contribution was given in the seventies by two groups which systemati- 
cally analysed the dispersion velocity in many spiral galaxies [3]: in fact, the velocity 
curves in the galaxy plane as a function of distance from the galactic center stay flat 
even outside the luminous disk, crediting the presence of a dark halo. Several other ex- 
perimental evidences for the Dark Universe have been pointed out by the progresses - 
with time passing - in the astronomical observations, such as: i) the Large Magellanic 
Cloud spins around our Galaxy faster than expected in case only luminous matter 
would be present; ii) the observation of X-ray emitting gases surrounding elliptical 
galaxies; iii) the velocity distribution of hot intcrgalactic plasma in clusters. All these 
observations have further supported that the mass of the Universe should be much 
larger than the luminous one in order to explain the observed gravitational effects. 

The existence of the Dark Universe is supported also by the standard cosmology 
(based on the assumption that the Universe arose from an initial singularity and went 
on expanding) in the inflationary scenario (proposed to avoid any fine tuning in the 
Big Bang initial conditions), which requires a flat Universe with density equal to the 
critical one: p c = |^ = 1.88ft- 2 • 10~ 29 g • cm~ 3 , where G is the Newton constant 
and H is the Hubble constant equal to 100ft. kms~ 1 Mpc _1 and 0.55 < h < 0.75. The 
uncertainty is due to the measurements of the actual value of the expansion rate of the 
Universe and to the considered models [4]; a recent determination from the WMAP 
data gives: ft = 0.72 ± 0.05 [5]. 

In particular, the density parameter fi = where p is the average density of 
the Universe (matter + energy), is a key parameter in the interpretation of the data 
from the measurements on Cosmic Microwave Background (CMB) since the global 
curvature of the Universe is related to it. The experimental results are consistent 
with a flat geometry of the Universe and, therefore, also support Q ~ 1 [6]; the most 
recent determination from the WMAP gives: il = 1.02 ± 0.02 [5]. Thus, the scenario 
is consistent with adiabatic inflationary models and with the presence of acoustic 
oscillations in the primeval plasma and requires the existence of Dark Matter in the 
Universe since the average density of the Universe as measured by photometric methods 
is: ft ~ 0.007. However, the detailed composition of ft in term of matter, O m , and of 
energy, Q^, cannot be inferred by CMB data alone; some information can be derived 
by introducing some other constraints [5, 7]. 

For the sake of completeness, we also mention that in last years studies have been 
performed [8] on astronomical standard candles as supernovae type la, that allow to 
evaluate relations between redshift and distance. These studies seem to point out an 
Universe whose expansion is accelerating, crediting the possible presence of a Dark 
Energy. When these results are combined with CMB data, f^A would account for 
about 70% of Q[5, 9]. This form of energy, with repulsive gravity and possible strong 
implication on the future evolution of Universe, would not be a replacement for Dark 
Matter and is still a mysterious task; dedicated ground and space based experiments 
are planned in order to confirm this scenario. 

Finally, as regards our Galaxy, from dynamical observations one can derive that it 
is wrapped in a dark halo, whose density nearby the Earth has been estimated to be 
for example in refs. [10, 11]: phaio — (0.17— 1.7) GcV cm~ 3 (see also later). 
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1.2 The nature of the Dark Matter 



The investigation on the nature of the Dark Universe has shown that large part of it 
should be in non-baryonic form. 

In fact, as regards baryons, in the past from the theory of big-bang nucleosynthesis 
(BBN) and from a lower limit to the primordial deuterium abundance a baryon density 
flB < 0.1 was set [12]. This upper limit has been precised by recent measurements 
of primordial deuterium abundance, giving f^ft 2 = 0.020 ± 0.001[13], that combined 
with the present determination of the Hubble constant implies: fts — 0.04; the latest 
determination by CBM experiments: Sl^/i 2 = 0.022 ± 0.003[5, 7], is also in good 
agreement. Recently, large efforts have been devoted to the investigation on Dark 
Baryonic Matter by experiments like EROS, MACHO and OGLE, which search for 
massive compact halo objects as baryonic candidates looking at microlensing effect 
toward Large and Small Magellanic Clouds and toward the Milky Way bulge. At 
present, in agreement with the expectations, the obtained results [14, 15] strongly 
limit the possible amount of Galactic Dark Matter in this form. In addition, a further 
argument, which also supports that the major part of the Dark Matter in the Universe 
should be in non-baryonic form, is the following: it is very difficult to build a model 
of galaxy formation without the inclusion of non-baryonic Dark Matter. 

Thus, a significant role should be played by non-baryonic relic particles from the 
Big Bang. They must be stable or with a lifetime comparable with the age of the Uni- 
verse to survive up to now in a significant amount. They must be neutral, undetectable 
by electromagnetic interactions and their cross section with ordinary matter should be 
weak (in fact, if their annihilation rate would be greater than the Universe expansion 
rate, they should disappear). The Dark Matter candidate particles are usually clas- 
sified in hot Dark matter (particles relativistic at decoupling time with masses < 30 
eV) and in cold Dark Matter (particles non relativistic at temperatures greater than 
10 4 K with masses from few GeV to the TeV region or axions generated by symmetry 
breaking during primordial Universe). The light neutrinos are the natural candidates 
for hot Dark Matter; they are strongly constrained by cosmology and a value over the 
limit f2„ ~ 0.05 gives an unacceptable lacking of small-scale structure[16]. In addition, 
a pure hot Dark Matter scenario is also ruled out by the measurements of the CMB 
radiation, which does not show sufficiently large inhomogeneity. 

Thus, cold Dark Matter candidates, which can be responsible for the initial grav- 
itational collapse, should be present and in large amount, although a pure cold Dark 
Matter scenario seems to be not favoured by the observed power spectrum of the den- 
sity perturbation. In practice, a mixed Dark Matter scenario is generally favourably 
considered. However, other possibilities can be considered such as, for example, the 
so-called "tilted Dark matter scenario" that introduces a significant deviation from 
the Zeldovich scale invariance of the power spectrum of the initial fluctuations. Any- 
how, in all the possible scenarios a significant fraction of cold Dark Matter particles is 
expected. 

As mentioned above, cold Dark Matter can be in form of axions or of WIMPs 
(Weakly Interacting Massive Particles). The axions are light bosons, hypothesized 
to solve the CP problem in strong interactions. Direct detection experiments are in 
progress since time by studying their interactions with strong electromagnetic fields, 
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Figure 1: Exclusion plot in the plane axion to photon coupling constant, ga-yy, versus 
axion mass, m a , achieved by DAMA/Nal in ref. [18]. The limit quoted in the paper 
(3077 < 1-7 x 10~ 9 GeU _1 at 90% C.L.) is shown together with the expectations of the 



KSVZ and DFSZ models; see ref. [18] for details. 



but no positive evidence has been found so far [17]. For completeness, we mention that 
some experiments (including DAMA/Nal, see Fig. 1) have also searched for possible 
axions produced in the Sun (see e.g. [18, 19]) and that some other will be realized in 
near future. However, these latter experiments cannot be classified as experiments for 
Dark Matter direct detection since they are not searching for relic axions. 

For the sake of completeness, we remind that also more exotic candidates (which 
generally could account for small fraction of Dark Matter in the galactic halo) have 
been considered and searched for, such as e.g. the magnetic monopoles with mass 
10 16 - 10 17 GeV [20], the neutral Strongly Interacting Massive particles (SIMPs) and 
the neutral nuclearities[21, 22, 23], the Q-balls[24], etc.; experimental searches for 
such candidates have given always negative results. Some of them have also been 
investigated by DAMA/Nal [22, 24]. 



2 The particles searched for 

The WIMPs are particles in thermal equilibrium in the early stages of the Universe, 
decoupled at freeze out temperature. Considering the WIMP particles as stable and 
with the same initial density for particles and antiparticles, their annihilation cross 
section, a an n, should be such that their annihilation rate should be lower than the 

expansion rate of the Universe: < a ann ■ v >~ ry-^ ^cmV 1 , where v is the 

relative velocity of the particle-antiparticle pair: thus, the interaction cross section 
is of the same order as those known of weak interactions. In case the particles and 
antiparticles would not have the same initial density, this relation would represent a 
lower limit. 



The velocity-spatial distribution of the WIMPs in our galactic halo is not well 
known. So far the simplest, non-consistent and approximate isothermal sphere model 
has generally been considered in direct WIMP searches; under this assumption the 
WIMPs form a dissipationless gas trapped in the gravitational field of our Galaxy in 
an equilibrium steady state and have a quasi-maxwellian velocity distribution with a 
cut-off at the escape velocity from the galactic gravitational field. More realistic halo 
models have been proposed by various authors such as Evans' power-law halos, Michic 
models with an asymmetric velocity distribution, Maxwellian halos with bulk rotation, 
etc. [25]. In particular, a devoted discussion on a wide (but still not complete) number 
of consistent halo models and their implications on available experimental data has 
been carried out e.g. in refs. [11, 25]; they will be summarized in §7.1.3. 

At present, the most widely considered candidate for WIMP is the lightest super- 
symmetric particle named neutralino, x- I n t ne Minimal Supersymmetric Standard 
Model (MSSM) where R-parity is conserved, the lightest SUSY particle, Xi must be 
stable and can interact neither by electromagnetic nor by strong interactions (other- 
wise it would condensate and would be detected in the galactic halo with the ordinary 
matter). The x is defined as the lowest-mass linear combination of photino (7), zino 
(Z) and higgsinos (hi, hi)- X = fl i7 + a iZ + a^h\ + 04/12 (where 7 and Z are lin- 
ear combination of U(l) and SU(2) neutral gauginos, B and W3) and is a Majorana 
particle. Under some assumptions, the x mass and the <n coefficients depend on the 
Higgs mass mixing parameter, fj,, on the B and W3 masses and on tg/3 (the ratio 
between the v.e.v's which give masses to up and down quarks). Thus, often the the- 
oretical estimates and sometimes the experimental results are presented in terms of 
(j,, tg(3 and wino mass, Mi. The \ cross section on ordinary matter is described by 
three Feynman diagrams: i) exchange between x an d quarks of the ordinary matter 
through Higgs particles (spin-independent - SI - interaction); ii) exchange between \ 
and quarks of the ordinary matter through Zo (spin-dependent - SD - interaction); 
iii) exchange between \ and quarks of the ordinary matter through squark (mixed - 
SI/SD - interaction). The evaluation of the expected rates for \ depends on several 
parameters and procedures, which are affected by significant uncertainties, such as 
e.g. the considered neutralino composition, the present uncertainties on the measured 
top quark mass and on certain sectors of the fundamental nuclear cross sections, on 
some lack of information about physical properties related to Higgs bosons and SUSY 
particles, on the possible use of constraints from GUT schemes and/or from b — > s + 7 
branching ratio, on the used rescaling procedure, etc.; in conclusion, considering also 
the large number of involved parameters, the supersymmetric theories have unlikely 
no practical predictive capability. 

Other candidates can also be considered as WIMPs; in particular, we remind an 
heavy neutrino of a 4-th family [26] and the sneutrino in the scenario described in ref. 
[27].' 

The heavy neutrino of a 4-th family was one of the first candidate proposed to 
solve the Dark Matter problem. Still now it may be considered as a good and realistic 
candidate, although unable to account for the whole Dark Matter missing mass. Such 
a neutrino could contribute - by its pair annihilation in the galactic halo - to positrons, 
antiprotons and diffused gamma background and these signatures might be better fit 
to the observed data [26]; moreover, it might dominate the Higgs decay mode in near 
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future LHC accelerator. The cosmological relic abundance of heavy neutrinos can be 
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Figure 2: Relic abundance of an heavy neutrino as a function of its mass according to 
the calculation of ref. [28]; masses above the Zo pole are considered. 

evaluated, as reported in Fig. 2, taking into account that the couplings are described 
within the Standard Model of elementary particles. Applying the condition that the 
density of such heavy neutrinos cannot exceed the critical density, a window in their 
mass can be evaluated [29]: 3 GeV < m v < 3 TeV. Considering the measurements 
of Z decay into invisible channels carried out at LEP and some implications of the 
measured cosmic ray flux [28], a mass range around 50 GeV with a reasonable local 
abundance (which permits to consider it as a Dark Matter candidate) is still open. 

In some supersymmetric models the lightest supersymmetric particle (LSP) can 
be the sneutrino, v, the spin-0 partner of the neutrino. In supersymmetric theories 
with no violation of leptonic number, a sneutrino with mass in the range 550 GeV < 

< 2300 GeV could have a relevant cosmological abundance (0.1 < Vl^h 2 < 1) [30]; 
however, because of its large interaction cross sections, the sneutrino cannot generally 
be considered as major component of Cold Dark Matter. Anyhow, a sneutrino as 
a candidate remains still possible in supersymmetric models with violation of lepton 
number [31]. In this framework the sneutrino can exist in two mass states, i>±, with 
a S ~ Am 2 /2m£, mass splitting (for Am 2 <C ml), being Am 2 a term introduced by 
the leptonic number violating operator. The two mass eigenstates have off-diagonal 
coupling with Zq boson and only couplings between i> + e i>- exist. As a consequence, 
the elastic scattering cross section on nuclei is extremely low [31] and sneutrinos with 
mass around 40-80 GeV and S about 5 GeV could have cosmological relic abundance 
in the range 0.1-1 [31]. Moreover, whatever scalars would be introduced in the the- 
ory, they can mix with sneutrinos and, consequently, the gauge interaction would be 
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reduced through the mixing angle [32] . The suppression of this interaction implies a 
sizeable relic abundance of the sneutrino even for low S values (e.g. around ~ 100 
keV). A similar sneutrino has been proposed as a possible WIMP candidate providing 
- through the transition from lower to upper mass cigcnstate - inelastic scattering 
with nuclei [27] (see also later). 

Finally, we remind that - in principle - even whatever massive and weakly inter- 
acting particle, not yet foreseen by theories, can be a good candidate as WIMP. 

In the following we will focus our attention on the WIMP direct detection technique 
in underground laboratory, where the low environmental background allows to reach 
the highest sensitivity; this is the process investigated by DAMA/Nal. We will later 
mention few arguments on the indirect detection approach, mainly in the light of some 
recent analyses. 

3 Some general arguments on the WIMP direct de- 
tection approach 

The WIMP direct detection approach mainly investigates the WIMP elastic scattering 
on the nuclei of a target-detector; the recoil energy is the measured quantity. In 
fact, the additional possibility to investigate the WIMP-nucleus inelastic scattering 
producing low-lying excited nuclear states (originating successive de-excitation gamma 
rays and, thus, presence of characteristic peaks in the measured energy spectrum) 
is disfavoured by the very small expected counting rate; for this reason, only few 
preliminary efforts have been carried out so far on this subject [33, 34, 35]. 

In the following subsections only few general arguments are addressed on the di- 
rect detection approach, while we simply remind that most experienced detection 
techniques have already been briefly commented in ref. [36], mainly in the light of a 
possible effective search for a WIMP signature. 

3.1 Some generalities 

A direct search for Dark Matter particles requires: i) a suitable deep underground site 
to reduce at most the background contribution from cosmic rays; ii) a suitable low 
background hard shield against electromagnetic and neutron background; iii) a deep 
selection of low background materials and a suitable identification of radio-purification 
techniques to build a low background set-up; iv) severe protocols and rules for building, 
transporting, handling, installing the detectors; v) an effective Radon removal system 
and control on the environment nearby the detectors; vi) a good model independent 
signature; vii) an effective monitoring of the running conditions at the level of accuracy 
required by the investigated WIMP signature. 

As an example of the suitable performances of a deep underground laboratory we 
remind those measured at the Gran Sasso National Laboratory of I. N.F.N, where the 
DAMA/Nal experiment has been carried out: i) muon flux: 0.6 muons m~ 2 h _1 [37]; 
ii) thermal neutron flux: 1.08 -10~ 6 neutrons cm~ 2 s _1 [38]; iii) epithermal neutron 
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flux: 1.98 -lCT 6 neutrons cm- 2 s" 1 ^]; iv) fast (E n > 2.5MeV) neutron flux: 0.09 
•10~ 6 neutrons cm" 2 s _1 [39]; v) Radon in the hall: ~ 10-30 Bq m" 3 [40]. 

The low background technique requires very long and accurate work for the se- 
lection of low radioactive materials by sample measurements with HP-Ge detectors 
(placed deep underground in suitable hard shields) and/or by mass spectrometer anal- 
yses; thus, these measurements are often difficult experiments themselves, depending 
on the required level of radiopurity. In addition, uncertainties due to the sampling pro- 
cedures and to the subsequent handling of the selected materials to build the apparata 
also require further time and efforts. As an example of an investigation of materials 
and detector radiopurity one can consider ref. [41], where the residual radioactivity 
measured in materials and detectors developed for DAMA/Nal is reported. More- 
over, some arguments on how to further improve the radiopurity of Nal(Tl) detectors 
(largely followed e.g. in the developments of the new DAMA/LIBRA set-up, now in 
test runs) can be found e.g. in ref. [42]. An interesting paper on the low background 
techniques is also e.g. ref. [43]. 

Main efforts regard the reduction of standard contaminants: 238 U and 232 Th (be- 
cause of their rich chains) and 40 K (because of its large presence in nature). When 
suitable radiopurity is reached for these components, the possible presence of non- 
standard contaminants should be also seriously investigated by devoted measurements. 
As shown e.g. in ref. [44] for the case of a ionizing Ge experiment, several orders of 
magnitude of rate reduction can be obtained with time and efforts in improving the 
experimental conditions. 

3.2 The "traditional" model dependent approach 

Since often the used statistics in direct experiments is very poor, the simple comparison 
of the measured energy distribution with an expectation from a given model framework 
is carried out. This "traditional" approach - the only one which can be pursued by 
cither small scale or very poor duty cycle experiments - allows only to calculate model 
dependent limits on WIMP-nucleus cross section at given C.L.. In fact, although for 
long time the limits achieved by this approach have been presented as robust reference 
points, similar results are quite uncertain not only because of possible underestimated 
systematics when relevant data handling and reduction is performed, but also because 
the result refers only to a specific model framework. In fact the model is identified 
not only by the general astrophysical, nuclear and particle physics assumptions, but 
also by the needed theoretical and experimental parameters and by the set of values 
chosen in the calculations for them. Some of these parameters, such as the WIMP 
local velocity, w , and other halo parameters, form factors' parameters, quenching fac- 
tor, etc. are also affected by significant uncertainties. Therefore the calculation of 
the expected differential rate, which has to be compared with the experimental one 
in order to evaluate an exclusion plot in the plane WIMP cross section versus WIMP 
mass, is strongly model dependent. As an example, Fig. 3 shows how an exclusion 
plot is modified by changing (within the intervals allowed by the present determina- 
tions) the values of the astrophysical velocities [45] . Analogous effects will be obtained 
when varying - within allowed values - every other of the several needed parameters 
as well as when varying every one of the general assumptions considered in the calcu- 
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Figure 3: Example of the effects due to the uncertainties in a given model framework 
when calculating exclusion plots. Here the simple case for the halo local velocity, 
v , and the escape velocity, v esc , is shown in case of spin-dependent coupled WIMPs 
as from ref. [45]. The top curve for each nucleus has been calculated - in a given 
model framework - assuming vo = 180 km/s and v esc — 500 km/s, while the lower 
one has been calculated assuming vo — 250 km/s and v esc = 1000 km/s; all the 
considered values are possible at present stage of knowledge. Analogous effects will 
be found for every kind of experimental result when varying experimental/theoretical 
parameters/assumption for whatever target-nucleus. 



lations. Thus, each exclusion plot should be considered only strictly correlated with 
the "cooking list" of the used experimental/theoretical assumptions and parameters as 
well as with detailed information on possible data reduction/selection, on efficiencies, 
calibration procedures, etc. Moreover, since WIMP-nucleus cross sections on different 
nuclei cannot directly be compared, generally cross sections normalized to the WIMP- 
nucleon one are presented; this adds further uncertainties in the results and in the 
comparisons, requiring the assumptions of scaling laws 3 . 

Thus, comparisons should be very cautious since they have not an universal char- 
acter. In addition, different experiments can have e.g. different sensitivity to the 
different possible WIMP couplings. 

In conclusion, this model dependent approach has no general meaning, no poten- 
tiality of discovery and - by its nature - can give only "negative" results. Therefore, 
experiments offering model independent signature for WIMP presence in the galactic 
halo are mandatory. 



3 We take this occasion also to stress that exclusion plots given in terms of cross sections on nucleus 
arc not model independent as quoted sometimes "traditionally" in literature, since they depend e.g. 
on the considered halo model, on the considered nuclear form factors, etc. 
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3.2.1 ... with electromagnetic background rejection technique 

In order to overcome the long and difficult work of developing very low background 
set-ups, strategies to reject electromagnetic background from the data are sometimes 
pursued. This can be realized in several scintillators by pulse shape discrimination 
(since electrons show a different decay time respect to nuclear recoils, as carried out 
in Nal(Tl) and LXe e.g. by DAMA/Nal in ref. [46] and by DAMA/LXe in ref. [47]) 
or by comparing, for the same event, two different signals (when the recoil/electron 
response ratio is expected to be different, such as heat/ionization in Ge or Si [48, 49] 
and heat /light in CaWC>4 [50, 51]). The first case offers a relatively safer approach 
than the second one since basic quantities (such as e.g. the sensitive volume) are 
well defined, while the second one is more uncertain. Just as an example, in case of 
heat/ionization read-out the precise knowledge of the effective sensitive volume for 
each one of the two signals and the related efficiencies as a function of the energy are 
required. A further discrimination strategy, which uses a two-phases gas/liquid Xenon 
detector with an applied electric field, has been also suggested for future experiments; 
there the light amplitudes of the primary and of the secondary scintillation pulses are 
compared [52]. However, in this case the discrimination critically depends e.g. on the 
definition of the real sensitive volume, on the dependence of the discrimination power 
with ionization position, on gas purity, etc. 

In every case, whatever strategy is followed, always only a statistical discrimination 
is possible (on the contrary of what is often claimed) because e.g. of tail effects 
from the two populations, from the noise, etc. Furthermore, the existence of known 
concurrent processes (due e.g. to end-range alphas, neutrons, fission fragments or in 
some case also the so-called surface electrons) , whose contribution cannot be estimated 
and subtracted in any reliable manner at the needed level of precision, excludes that an 
unambiguous result on WIMP presence can be obtained following a similar approach. 

Moreover, when using similar procedures, the real reached sensitivity is based e.g. 
on the proper estimate of the systematic errors, on the accuracy of all the involved 
procedures and on the proper accounting of all the related efficiencies, on the proper 
knowledge of the energy scale and energy threshold (see also §7.1.6) and on the verified 
stability of the running conditions. Consider e.g. the difficulty to manage the efficiency 
due to the coincidence of the few keV heat/ionization or heat/scintillation signals or, 
in case of the two-phases LXe detectors, the triggering of the primary and secondary 
scintillations. We note also that sometimes in literature some methodologically uncor- 
rect methods are also considered which allow to claim for a larger sensitivity than the 
correct one. 

In conclusion, the possibility to achieve a control of the systematic error in rejection 
procedures at level of ~ 10~ 4 , as it has recently been claimed (see §5.1.1), appears 
unlikely whatever rejection approach would be considered. 

Finally, it is worth to note that rejection strategies cannot safely be applied to the 
data when a model independent signature based on the correlation of the measured 
experimental rate with the Earth galactic motion is pursued (see later); in fact, the 
effect searched for (which is typically at level of few %) would be largely affected by 
the uncertainties associated to the - always statistical - rejection procedure. On the 
other hand the signature itself acts as an effective background rejection as pointed out 
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e.g. for the WIMP annual modulation signature since ref. [53]. 

3.3 An unambiguous signature for WIMPs in the galactic halo 
is needed 

To obtain a reliable signature for WIMPs is necessary to follow a suitable model inde- 
pendent approach. In principle, three main possibilities exist; they are based on the 
correlation between the distribution of the events, detected in a suitable underground 
set-up, with the galactic motion of the Earth. 




Figure 4: Left: schematic representation of the experimental approach considered in 
ref. [54] to investigate the correlation between the recoil direction and the Earth ve- 
locity direction by using anisotropic scintillators. The anisotropic scintillator is placed 
ideally at LNGS with d axis in the vertical direction and b axis pointing to the North. 
The area in the sky from which the WIMPs are preferentially expected is highlighted. 
Right: expected rate, in the 3-4 keV energy window, versus the detector (or Earth) 
possible velocity directions. This example refers to the particular assumptions of a 
WIMP mass equal to 50 GeV, a WIMP-proton cross section equal to 3 • 10~ 6 pb and 
to the model framework of ref. [54]. The dependence on the "polar-azimuth" angle 
(<j> P a) induces a diurnal variation of the rate. 

The first one correlates the recoil direction with that of the Earth velocity, but 
it is practically discarded mainly because of the technical difficulties in reliably and 
efficiently detecting the short recoil track. Few R&D attempts have been carried 
out so far such as e.g. [55, 56], while a suggestion - based on the use of anisotropic 
scintillators - was originally proposed by DAMA collaborators in ref. [57] and recently 
revisited in ref. [54]. As an example, Fig. 4 (left) shows a schematic representation of 
the experimental approach studied in ref. [54] ; an example of the dependence of the 
expected rate on the WIMP arrival direction, with respect to the crystal axes, for the 
considered experimental case is given in Fig. 4 (right). 
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The second approach correlates the time occurrence of each event with the diurnal 
rotation of the Earth. In fact, a diurnal variation of the low energy rate in WIMP 
direct searches can be expected during the sidereal day since the Earth shields a 
given detector with a variable thickness, eclipsing the WIMP "wind" [58]. However, 
this effect can be appreciable only for relatively high cross section candidates and, 
therefore, it can only test a limited range of Cold Dark Matter halo density. For a 
recent experimental result see e.g. ref. [59], where a statistics of 14962 kg-day collected 
by DAMA/Nal has been investigated in the light of this signature. As an example the 
dependence of 9 (the angle defined by the Earth velocity in the Galactic frame with 
the vector joining the center of the Earth to the position of the laboratory) on the 
sidereal time, is shown in Fig. 5 (left) in case of the Gran Sasso National Laboratory 
location. The expected signal rate, in case of the experimental set-up and assumptions 
quoted in ref. [59], is given in Fig. 5 (right). 




Figure 5: Schematic description of the approach which correlates the time occurrence 
of each event with the diurnal rotation of the Earth. Left: the 9 angle (defined by the 
Earth velocity in the Galactic frame with the vector joining the center of the Earth 
to the position of the laboratory) as a function of the sidereal time; here the case 
for the Gran Sasso National Laboratory of the I. N.F.N, is considered. Right: signal 
rate expected in the 2-6 keV energy interval when assuming a 60 GeV WIMP mass, 
a WIMP-proton cross section equal to: a) 7.0 • 10~ 6 pb, b) 5 • 10~ 2 pb, c) 10 _1 pb, d) 
1.0 pb, and the model framework of ref. [59]. 

The third possibility, feasible and able to test a large interval of cross sections and 
of WIMP halo densities, is the so-called annual modulation signature [53]. This is the 
main signature exploited by DAMA/Nal [60, 61, 62, 63, 64, 65, 66, 11]. The annual 
modulation of the signal rate is induced by the Earth revolution around the Sun; as a 
consequence, the Earth is crossed by a larger WIMP flux in June (when its rotational 
velocity is summed to the one of the solar system with respect to the Galaxy) and by 
a smaller one in December (when the two velocities are subtracted) (see Fig. 6). 

In particular, the expected differential rate as a function of the recoil energy, 
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-30 km/s 



Figure 6: Schematic view of the Earth motion around the Sun. 



dR/dE R (see §7.1 for detailed discussion), depends on the WIMP velocity distribution 
and on the Earth's velocity in the galactic frame, v e (t). Projecting v e (t) on the galactic 
plane, one can write: 

V e (t) — Vq + V^COS^COSUj{t — to) (1) 

here Vq is the Sun's velocity with respect to the galactic halo (u© ~ vq + 12 km/s and 
v is the local velocity whose value is in the range 170-270 km/s [62, 67]); v© = 30 
km/s is the Earth's orbital velocity around the Sun on a plane with inclination 7 = 
60° respect to the galactic plane; furthermore, u>= 2tt/T with T=l year and roughly 
to ~ 2 nd June (when the Earth's speed is at maximum). The Earth's velocity can be 
conveniently expressed in unit of "Uo: i](t) — v e (t)/vo = rjo + Arjcosuj{t — to), where 
- depending on the assumed value of the local velocity - ?7o=1.04-1.07 is the yearly 
average of 77 and Ar] = 0.05-0.09. Since A77 <C r)o, the expected counting rate can be 
expressed by the first order Taylor approximation: 



dR , , „ dR , , d 



I''!"] = TEH'/"] + — \4Sz J ArjcosLuit-to). (2) 

no 



Averaging this expression in a A:-th energy interval one obtains: 

Sk[w(t)] = Skblo] + [-^—\vo^'n cos<jJ {t - t ) = S ,k + S m ,kCosu(t - to), (3) 

with the contribution from the highest order terms less than 0.1%. The first time- 
independent term is: 

So,k = ^r[ -^[vo]dE R , (4) 
AE k JAE k aE R 

while the second term is the modulation amplitude given by: 

max mini / r \ 
S m .k = TTT / TT- "Trr! Ar]dE R ~ , 5 

A^fe J AEk drj \dE R J 71=Vo 2 
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with rjmax = r] + At] and rj min = rjo - Arj. The 5o,fc and S m ,k are functions of 
the parameters associated with the WIMP interacting particle (such as e.g. mass 
and interaction cross sections), of the experimental response of the detector, of the 
considered model framework and of the related parameters (see later). 

It is worth to note that the S mi £ values can be not only positive, but also negative 
or zero, due to the expected energy distribution profiles in June and in December 
within a finite energy window [68]. Therefore, the highest sensitivity can be obtained 
when considering the smallest energy bins allowed by the available statistics in the 
energy region of interest. 

Although the modulation effect is expected to be relatively small (the fractional 
difference between the maximum and the minimum of the rate is of order of ~ 7%), 
a suitable large-mass, low-radioactive set-up with an efficient control of the running 
conditions - such as DAMA/Nal [41] - would point out its presence. In fact, a suit- 
able correlation analysis can allow to extract even a small periodic component, su- 
perimposed with a time independent signal and a background [53] . With the present 
technology, the annual modulation remains the main signature of a WIMP signal. 

In addition, the annual modulation signature is very distinctive since a WIMP- 
induced seasonal effect must simultaneously satisfy all the following requirements: the 
rate must contain a component modulated according to a cosine function (1) with one 
year period (2) and a phase that peaks roughly around ~ 2 nd June (3); this modulation 
must only be found in a well-defined low energy range, where WIMP induced recoils 
can be present (4); it must apply to those events in which just one detector of many 
actually "fires", since the WIMP multi-scattering probability is negligible (5); the 
modulation amplitude in the region of maximal sensitivity must be <7% (6). Only 
systematic effects able to fulfil these 6 requirements could mimic this signature and 
- as far as we know - no other effect investigated so far in the field of rare processes 
offers a so stringent and unambiguous signature. 

Of course, the amount of the measured effect depends e.g. on the sensitivity of the 
experiment to the coupling of the WIMP candidate, on the WIMP particle physics 
features, on the nuclear features of the used target- nucleus and on the quality of the 
running conditions. 

4 The DAMA experiment 

The DAMA experiment has been worked and works as an observatory for rare processes 
(such as WIMP direct detection, (3(3 decay processes, charge-non-conserving processes, 
Pauli exclusion principle violating processes, nucleon instability, solar axions and ex- 
otics [18, 22, 24, 34, 35, 41, 46, 47, 59, 60, 61, 62, 63, 64, 65, 66, 68, 69, 70, 71, 72, 11]) 
by developing and using low radioactive scintillators. It is installed deep underground 
in the Gran Sasso National Laboratory of I.N.F.N.. 

The main developed and used experimental set-ups are: the ~ 100 kg Nal(Tl) 
set-up (DAMA/Nal) [41] (which has completed its data taking in July 2002), the 
~ 6.5 kg liquid Xenon set-up (DAMA/LXe) [70], the so-called "R&D" apparatus 
(DAMA/R&D) and the new LIBRA (Large sodium Iodide Bulk for RArc processes; 
~ 250 kg of ultra-radiopure Nal(Tl)) set-up (DAMA/LIBRA) whose installation has 
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been started in fall 2002 and which is presently in test run. Moreover, an underground 
low-background germanium detector allows to select materials for radiopurity. 

In the following the final model independent result of DAMA/Nal on the inves- 
tigation of the WIMP annual modulation signature is discussed. We remind that 
DAMA/Nal is the largest mass, highest sensitivity experiment, built before the new 
DAMA/LIBRA, having as main aim the investigation of WIMPs in the galactic halo. 
It was a pioneer experiment proposed in 1990 [73], which has opened for other experi- 
ments and approaches in the field; moreover, its results on the investigation of WIMPs 
in the galactic halo by the annual modulation signature have - by the fact - motivated 
the wide interest in the field arisen in recent years. 

4.1 DAMA/Nal 

The DAMA/Nal set-up [41] can effectively exploit the WIMP annual modulation sig- 
nature because of its well known technology, of its high intrinsic radiopurity, of its 
mass, of its suitable control of all the operational parameters and of the deep under- 
ground experimental site. 

The detailed description of the ~ 100 kg Nal(Tl) DAMA set-up, of its radiopu- 
rity, of its performance, of the used hardware procedures, of the determination of the 
experimental quantities and of the data reduction has been given in refs. [41, 63, 64]. 

Here we only recall that the detectors used in the annual modulation studies are 
nine 9.70 kg highly radiopurc Nal(Tl) scintillators especially built for this purpose 
in a joint effort with Crismatec company. The bare Nal(Tl) crystals are encapsu- 
lated in suitably radiopure Cu housings; 10 cm long Tetrasil-B light guides act as 
optical windows on the two end faces of the crystals and are coupled to specially 
developed EMI9265-B53/FL photomultipliers (PMT), which are supplied by positive 
voltage with ground cathode. The two PMTs of a detector work in coincidence and 
their threshold is set at the single photoelectron level; the measured light response 
is 5.5 - 7.5 photoelectrons/keV depending on the detector [41]. The software energy 
threshold has been cautiously taken at 2 keV [41, 46, 64]. The detectors are inside 
a low radioactivity sealed copper box installed in the center of a low radioactivity 
Cu/Pb/Cd-foils/polyethylene/paraffin shield. Moreover, about 1 m concrete (made 
from the Gran Sasso rock material) almost fully surrounds (outside the barrack) this 
passive shield, acting as a further neutron moderator. The copper box is maintained 
in a high purity (HP) Nitrogen atmosphere in slightly overpressure with respect to the 
external environment. Furthermore, also the whole shield is sealed and maintained in 
the HP Nitrogen atmosphere. The whole installation is air-conditioned and the tem- 
perature is monitored. On the top of the shield a glove-box (also maintained in the HP 
Nitrogen atmosphere) is directly connected to the inner Cu box, housing the detectors, 
through Cu pipes. The pipes are filled with low radioactivity Cu bars (covered by 10 
cm of low radioactive Cu and 15 cm of low radioactive Pb) which can be removed 
to allow the insertion of radioactive sources for calibrating the detectors in the same 
running condition, without any contact with external air [41]. 

An hardware/software system to monitor the running conditions has been opera- 
tive; in particular, several probes have been read out by the data acquisition system 
and stored with the production data. Moreover, self-controlled computer processes are 
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operational to automatically control several parameters and to manage alarms [41, 64]. 

The electronic chain and the data acquisition system used during the DAMA/Nal- 
to DAMA/NaI-5 running periods has been described in ref. [41]. At completion of 
the DAMA/NaI-5 data taking (summer 2000) the whole electronics and DAQ have 
been completely substituted; they are briefly summarized in the following. This new 
system has been operative during DAMA/NaI-6 to -7 running periods, that is up to 
the end of the DAMA/Nal data taking. 

The new DAQ system has been based on a Digital Alpha Workstation with Dig- 
ital Unix operating system interfaced with the VXI and CAMAC components of the 
electronic chain via a GPIB bus; the acquisition program has been developed on the 
basis on the system discussed in ref. [74] and on a specific applicative software. In 
the new configuration the HV power supply for the PMTs has been given by a CAEN 
multichannel voltage supply with voltage stability of 0.1%. 



Preamp 



Main Trigger 


L19 


Gate I 

ADC 1 L 


gate 1 






200re 


CAMAC r 


J£. LB 





pfflf' 



Preamp 
Gain 10 





WA1 VXI 


f[fT] L16 ' 


In 1 



OR enable 1 
Line CAMAC 



Scaler 1 
CAMAC 



L24 J 

L25 


Single Trigger 
crystal 1 


L2G 

L27 


Single Trigger 
crystal 1 



| L23ybate 1 f | 


IR pattern 








L28 


recognition 




In 1 tJ 





Delay line 
200rs 


ADC 2 L 
CAMAC r 

Gate 


Main Trigger 
gate 2 


L20 








pattern ~~j 




alog 5 frretal 9 1 — ^ I 



Figure 7: The new electronic chain installed after the completion of the DAMA/Nal- 
5 running period. Top: schema of the electronic chain of a single detector with its 
trigger in the new electronic chain. Bottom left: schema of the main trigger of the 
new acquisition system. Bottom right: schema of the trigger of the new Waveform 
Analyzers (WM.). 
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In Fig. 7 the analogic part of this new electronic chain for one detector and its 
trigger are shown as well as the main trigger of the acquisition system and the trigger 
system of the new waveform analyzer. We note that the analogic part and the trigger 
for single detector are similar to the ones in the previous electronic chain, but the 
multiplexer system - which was used in the past (having there at disposal only one 
single channel Transient Digitizer LeCroy 8828D with 200 MSamplc/s sampling fre- 
quency) - has been removed since now each detector has a devoted waveform analyzer 
(WA) channel. This is accomplished using fast VXI Tektronix four-channel TVS641A 
digitizers with a sampling frequency of 1 GSample/s and 250 MHz bandwidth. The 
digitizers provide a VXI word-serial protocol for commands, while a specific Fast Data 
Channel (FDC) driver developed by [75] has been used to speed up the data transfer 
through the GPIB bus. The main trigger part of the electronic schema and the high 
energy data acquisition are also similar to the ones in the previously used electronic 
chain [41]. 

For completeness and for template purpose, let us briefly describe this new elec- 
tronic chain. The signal LI (see Fig. 7 for the definition of the symbols) from the first 
PMT is sent to the Al preamplifier having 0-250 MHz bandwidth, a factor 10 gain 
and a voltage integral linearity ±0.2%. The signal L2 from the second PMT is divided 
in two branches: 19/20 of the signal is sent to the input of the A2 preamplifier, while 
the remaining 1/20 - suitably delayed - feeds a charge ADC (ADC3) channel. This 
last part processes the pulses with amplitude such to saturate the remaining part of 
the electronics (they correspond to high energy events). The preamplified signals - 
through linear Fan-in/Fan-out devices - provide the inputs for the charge ADCs (left 
signals and right signals) and for the VXI waveform analyzer (WM1 in Fig. 7) (which 
processes the signal in a 2 us time window) 4 . 

The electronic devices, that provide the trigger of a single detector, are shown in 
Fig. 7. In particular, the copies of the PMT signals arc the inputs of the Timing Filter 
Amplifiers (TFA1 and TFA2) which amplify and integrate the signal (integration 
time 50 ns); their outputs are discriminated (Discriml and Discrim2) with single 
photoelectron trigger level. The coincidence (AND1) between the two logical NIM 
outputs provides the single trigger of the detector. The other inputs of AND1 are: 
the signal L23 given by a Gate Generator (GG1) which allows to reject afterglow and 
Bi-Po events in a 500/xs time window after the occurrence of the event (introducing 
a systematic error on the measured rate of about 10~ 4 ) and the signal L22 given by 
a CAMAC I/O Register. The latter permits to enable or disable the single detector 
trigger during the calibrations. The outputs of ANDl provide: i) the signal for a 
CAMAC scaler (Scalerl) to count the events for each detector; ii) the L25 and L26 
used in the main trigger (see later); iii) L27 giving the start to the Gate Generator 
(GG1) which - in addition to the veto of the coincidence - gives the signal L28 issued 
to a 16-bit CAMAC I/R Pattern Recognition which allows to identify the detector or 
the detectors which have generated the trigger. 

The general trigger of the acquisition - see Fig. 7 - is provided by the logic OR1 of 
all the crystals. The output of AND14 is issued to a Scaler, which counts the number 

4 For completeness we note that the 4 detectors named SIMP in rcf. [41] which have been used 
only as additional anticoincidence and - sometimes - in special triggers, have not been equipped in 
the new electronic chain with waveform analyzers. 
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of main triggers sent to the coincidence AND15. The coincidence AND15 generates 
triggers only when the acquisition is ready to manage them. Therefore, the dead time 
of the acquisition is properly accounted in the estimate of the running time by using 
the information from the scaler after AND14. When a general trigger occurs, the 
following logic signals are issued to: i) the Gate Generator (GG2) generating the 600 
ns gates the charge ADCs; ii) the Delay Gate Generator which gives the strobe signal 
to the I/R Pattern Recognition and generates the LAM (and, therefore, the interrupt 
to the CPU of the acquisition computer) in the CAMAC system; iii) the Delay Gate 
Generator which gives the signal to the trigger of the waveform analyzers. This last 
condition is verified only if the total energy deposited in the detectors is in an energy 
window suitably chosen (1 to 90 keV). For this purpose, each line feds a Spectroscopy 
Amplifier whose gain is equalized in order to have the same response for each detector. 
Therefore, a Single Channel Analyzer made by the two discriminators, Discrim LL 
and Discrim UL, allows to select only events in the chosen energy window. 

A devoted electronic circuit [76], shown in Fig. 7, allows to trigger only the Wi's 
which correspond to fired detectors; it gives a trigger to each WA when: i) at least 
one of its corresponding lines has a trigger; ii) the main trigger is present; iii) the total 
energy of the events is in the chosen energy window. Let us remind that for the events 
with energy outside this energy window (e.g. high energy events) the ADC values are 
acquired in any case. 

As regards other aspects, we recall that the linearity and the energy resolution of 
the detectors have been investigated using several sources [41, 64] such as, for the low 
energy region, 55 Fe (5.9 keV X-rays), 109 Cd (22 keV X-rays and 88 keV 7 line) and 
241 Am (59.5 keV 7 line) sources. In particular, in the production runs, the knowledge 
of the energy scale is assured by periodical calibrations with 241 Am source and by 
monitoring (in the production data themselves summed every ~ 7 days) the position 
and energy resolution of the 46.5 keV 7 line of the 210 Pb [41, 46, 60, 61, 63, 64]. The 
latter peak is present - at level of few counts per day per kg (cpd/kg) - in the measured 
energy distributions mainly because of a contamination (by environmental Radon) of 
the external surface of the crystals' Cu housings, occurred during the first period of the 
underground storage of the detectors. The calibration sources are introduced in the 
proximity of the detectors by means of the pipes connected - as already described - 
with the upper " glove-box" , which is also continuously maintained in the HP Nitrogen 
atmosphere. 

As in every experiment in the field, obvious noise events (whose number sharply 
decreases when increasing the number of available photoclcctrons) have to be removed; 
the used procedure has been described e.g. in refs. [41, 64] 5 . We remind that the 
noise in this experiment is given by PMT fast single photoclcctrons with decay times 
of the order of tens ns, while the "physical" (scintillation) pulses have decay times 
of order of hundreds ns. Thus, the large difference in decay times and the relatively 
large number of available photoelectrons response assure an effective noise rejection 

5 This procedure assures also the rejection of any possible contribution cither from afterglows (when 
not already excluded by the dedicated 500 [is veto time; see above) induced by high energy events 
or from any possible Cerenkov pulse in the light guide or in the PMTs; in fact, they also have time 
decay of order of tens ns as the noise events. 
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[41, 63, 64]. Several variables can be built by using the pulse information recorded 
by the waveform analyzer [41, 64]). In particular, for each energy bin, we plot the 

V _ Areajfrom ns to 50 ns) i ....... t l.„ y _ Areajfrom 100 ns to 600 ns) i 

1 - Area(from ns to 100 ns) VaiUe verSUS tne A ~ Areatfrom ns to 600 ns) VaiUe 

calculated for every event. In the X, Y plane the slow scintillation pulses are grouped 
roughly around (X ~ 0.7, V ~ 0.5) well separated from the noise population which is 
grouped around small X and high Y values (sec e.g. rcf. [41]). The scintillation pulses 
are selected by applying an acceptance window in X, Y. Since the statistical spread of 
the two populations in the X, Y plane becomes larger when the number of available 
photoclcctrons and the signal/noise ratio decrease, windows with smaller acceptance 
become necessary to maintain the same noise rejection power. In the DAMA/Nal 
experiment they are kept enough stringent to assure also the absence of any possible 
residual noise tail in the scintillation data to be analysed [41]. According to standard 
procedures, the acceptance of the considered window for scintillation pulses in the A, Y 
plane is determined by applying the same procedure to the scintillation data induced 
- in the same energy intervals - by calibration sources [41, 63, 64]. In particular, for 
this purpose, about 10 4 - 10 5 events per keV are typically collected in the low energy 
region just above the 2 keV software energy threshold during routine calibration runs 
[41, 64]. All the periodical long calibration procedures [41, 46] and the time specifically 
allocated for maintenance and/or for improvements are the main components affecting 
the duty cycle of the experiment. Moreover, in the DAMA/NaI-1 running period the 
data have been taken only in the two extreme conditions for the annual modulation 
signature (see Table 1). 

The energy threshold, the PMT gain, the electronic line stability are continuously 
verified and monitored during the data taking by the routine calibrations, by the 
position and energy resolution of the 210 Pb line (see above) and by the study of the 
hardware rate behaviours with time. 

In particular, the measured low energy distributions of interest for the WIMP 
investigation have been given in refs. [59, 63, 64, 71], where the corrections for effi- 
ciencies and acquisition dead time have already been applied. We note that usually in 
DAMA/Nal the low energy distributions refer to those events where only one detector 
of many actually fires (that is, each detector has all the others in the same installation 
as veto; this assures a background reduction, which is of course impossible when a 
single detector is used). 

5 The first DAMA/Nal results on the annual mod- 
ulation signature 

The presence of a model independent effect has been firstly pointed out since the 
TAUP conference in 1997 [77] and corollary model dependent quests for a candidate 
particle have been analysed in some of the many possible model frameworks, improv- 
ing the quest with time (see Table 1). Cumulatively during four annual cycles a model 
independent effect (exposures up to 57986 kg • day; see Table 1) has been pointed 
out [63, 64]. No systematics or side reactions able to mimic the annual modulation 
signature has been found [64] ; this can be well understood when considering the partic- 
ularly stringent and numerous specific requirements which identify the WIMP annual 
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modulation signature itself (sec §3.3). No other experiment has at present suitable 

Table 1: Summary of the first running periods which have already cumulatively shown 
a 4 (j C.L. model independent effect. The related references are given in the third 
column. In the last column the improvements with time in the model-dependent quest 
for the candidate are summarized. 



Periods 


Statistics (kg x day) 


Rcf. 


Considered scenarios 
in corollary quests for the candidate 
(prior on mw from 
accelerator expts included) 


DAMA/NaI-1 


3363.8 (winter) 
+ 1185.2 (summer) 


[60] 


WIMP pure SI, 
Isothermal spherical halo, 
v = 220 km/s, 
Helm Form Factor, 
All the parameters fixed 
to their central value 


DAMA/NaI-2 


14962 
(Nov. — > end of July) 


[61], 
[62] 


as DAMA/NaI-1 
+ halo (co-)rotation 
+ uncertainty on vq 


DAMA/NaI-3 


22455 

(middle Aug. — * end of Sept.) 


[63] 


as DAMA/NaI-2 
+ prior from DAMA/NaI-0 


DAMA/NaI-4 


16020 

(middle Oct. — > middle Aug.) 


[63] 


as DAMA/NaI-3 


TOTAL 


57986 


[63], 
[64], 
[65], 
[66], 
[11] 


as DAMA/NaI-3 and DAMA/NaI-4 
+ SI&SD, "inelastic" 
+ Other consistent halo models, 

SD form factor from [117] 
Uncertainties on some parameters 


+ DAMA/NaI-0 


4123.2 


[46] 


Limits on recoils measured 
by pulse shape 
discrimination 



sensitivity, mass and control of the running conditions to effectively exploit the WIMP 
annual modulation signature as DAMA/Nal. 

As mentioned, the implications of the observed model independent effect have been 
in addition studied also under some - of the many possible - different model-dependent 
frameworks. In fact, some scenarios for purely spin-independent (SI), purely spin- 
dependent (SD) 6 , mixed SI and SD coupled WIMPs and also WIMPs with preferred 
inelastic scattering [60, 61, 62, 63, 64, 65, 66, 11] have been considered, including in 
the data analyses the lower bound on the mass of the supersymmetric candidate as 
derived from the LEP data in the usually adopted supersymmetric schemes based on 
GUT assumptions as e.g. in ref. [78]. This corollary investigation on the quest for 
a candidate particle has been improved with time in several aspects as summarized 
in Table 1. Theoretical implications of these results in terms of a neutralino with 
dominant SI interaction have been discussed e.g. in ref. [79, 80] for some theoretical 
model frameworks and in terms of an heavy neutrino of the fourth family in ref. [26] . 



6 For the sake of completeness, we comment that JHEP 0107 (2001) 044 is not at all in conflict 
with a possible SD solution since it considered only two particular purely SD couplings (of the many 
possible) in a strongly model dependent context and using modulation amplitudes valid instead only 
in a particular purely SI case. 
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5.1 Comparison with some model dependent results 
5.1.1 ... from direct searches 

As mentioned above no other experiment directly comparable with the model inde- 
pendent DAMA/Nal result on WIMPs in the galactic halo is available at present. 

Only few experiments [48, 49, 81], which use different target nuclei and different 
methodological approaches, have released extremely poor statistics following the so- 
called model dependent "traditional" approach (see section 3.2). We have reported 
several times (see e.g. rcf. [82]) some specific arguments; here we only summarize in 
Table 2 some main items. 

In particular, these experiments exploit a huge data selection releasing typically ex- 
tremely poor exposures with respect to generally long data taking and, in some cases, 
to several used detectors. Their counting rate is very high and few/zero events are 
claimed after applying several strong and hardly safe rejection procedures (involving 
several orders of magnitude; see Table 2). These rejection procedures are also poorly 
described and, often, not completely quantified. Moreover, most efficiencies and phys- 
ical quantities entering in the interpretation of the claimed selected events (see Table 
2) have never been discussed in the needed details; as an example, we mention the 
case of the quenching factor of the recoil target nuclei in the whole bulk material for 
the bolometer cases, which is arbitrarily assumed to be 1 (see §7.1.6), implying a sub- 
stantially arbitrarily assumed energy scale and energy threshold. The reproducibility 
of the results over different running periods has also not been proved as well as the 
values of the effective sensitive volumes for read-outs of the two signals (when applied) 
and the overall efficiencies. Further uncertainties are present when, as in ref. [48], a 
neutron background modeling and subtraction is pursued in addition. 

As regards in particular the Zeplin-I result of ref. [81, 83], a very low energy 
threshold is claimed (2 keV), although the light response is very poor: between ~ 1 
ph.e./keV [81] (for most of the time) and ~ 2.5 ph.e./keV (claimed for 16 days) [83] 7 . 
Moreover, a strong data filtering is applied to the high level of measured counting rate 
(~ 100 cpd/kg/keV at low energy, which is nearly two orders of magnitude larger that 
the DAMA Nal(Tl) background in the same energy region) by hardware vetoes, by 
fiducial volume cuts and, largely, by applying down to few keV a standard pulse shape 
discrimination procedure, although the LXe scintillation pulse profiles (pulse decay 
time < 30 ns) are quite similar even to noise events in the lower energy bins and in 
spite of the poor light response. Quantitative information on experimental quantities 
related to the used procedures has not yet been given [81, 83] 

In addition to the experimental aspects, these experiments generally perform an 
uncorrect quotation of the DAMA/Nal first quests for a purely SI coupled candi- 
date in some given model frameworks and ignore the published interpretation of the 
DAMA/Nal model independent effect in terms of candidates with other kind of cou- 
plings. Anyhow, intrinsically no reliable result can be achieved in a comparison of the 
exclusion plots quoted by these experiments with regions allowed by DAMA/Nal in 
corollary quests for a candidate. In fact, any exclusion plot always refers to a partic- 

7 For comparison we remind that the data of the DAMA/LXe set-up, which has a similar light 
response, are analysed by using the much more realistic and safer software energy threshold of 13 keV 
[47]. 
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Table 2: Features of the first DAMA/Nal results on the WIMP annual modulation 
signature (57986 kg x day exposure) [60, 61, 62, 63, 64, 65, 66, 11] with those of refs. 
[48, 49, 81]. See text. Here (as well as in the text) keV always means keV electron 
equivalent if not otherwise mentioned. 





DAMA/Nal 


CDMS-I 


Edclwciss-I 


Zeplin-I 


Signature 


annual 
modulation 


None 


None 


None 


Target-nuclei 


23 Na, 127 I 


"«'Ge 




nat Xe 


Technique 


well known 


poorly 
experienced 


poorly 
experienced 


critical optical 
liquid/gas interface 
in this realization 


Target mass 


~ 100 kg 


0.5 kg 


0.32 kg 


~ 3 kg 


Exposure 


57986 kg x day 


15.8 kg x day 


8.2 kg x day 


280 kg X day 


Depth of the 
experimental site 


1400 m 


10 m 


1700 m 


1100 m 


Software energy- 
threshold 


2 keV 
(5.5 - 7.5 p.c./kcV) 


10 keV 


20 keV 


2 keV 
(but: a/E = 100% 
mostly 
1 p.e./keV; [81]) 

(2.5 p.e./kcV 
for 16 days; [83]) 


Quenching 
factor 


Measured 


Assumed — 1 


Assumed — 1 


Measured 


Measured event 
rate in low 
energy range 


~ 1 cpd/kg/kcV 


~ 60 cpd/kg/kcV 
(10 5 events) 


2500 events 
total 


~ 100 cpd/kg/keV 


Claimed events 
after rejection 
procedures 




23 in Gc, 4 in Si, 
4 multiple evts in Ge 
+ MontcCarlo on 
neutron flux 





~ 20-50 cpd/kg/keV 
after rejection and 
?? after standard PSD 
[81, 83] 


Events satisfying 
the signature 
in DAMA/Nal 


modulation 
amplitude 
integrated over the 
given exposure 
2000 events 








Expected number 

of events from 
DAMA/Nal effect 




from few down 
to zero depending 

on the models 
(and on quenching 
factor) 


from few down 
to zero depending 

on the models 
(and on quenching 
factor) 


depends on 
the models 
(even zero) 



ular model framework where, in addition, all the involved nuclear cross sections are 
scaled to cross section(s) on nucleon (see §7.1); thus, it has no "universal" validity and 
- even within the same general assumptions for a model (e.g. purely SI coupling) - 
the proper accounting for parameters uncertainties, scaling laws uncertainties, form 
factors uncertainties, halo model uncertainties, etc. (see §7.1.3) significantly vary the 
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result of any comparison (even when assuming as correct the evaluation of the selected 
number of events and the energy scale and energy threshold determinations given in 
refs. [48, 49, 81, 83]). Moreover, there exist scenarios (see e.g. later in §7.1.3) to which 
Na and I are sensitive and other nuclei, such as e.g. no *Ge, ™ at Si and ™ at Xe, are not. 
Just as an example, a possible WIMP with a SI cross section of few 10~ 7 pb and with 
SD cross section of few 10 _1 pb would produce a sizeable signal in DAMA/Nal, but 
almost nothing in the Ge and Si experiments of refs. [48, 49] as well as in the Xenon 
target of ref. [81, 83] if the SD component would have ~ or 9 ~ 7r (see §7.1.1). 
In conclusion: 

1. no other experiment, whose result can be directly comparable in a model inde- 
pendent way with that of DAMA/Nal, is available so far. 

2. as regards in particular CDMS-I, EDELWEISS-I and Zeplin-I, e.g.: 

i) they are insensitive to the model independent WIMP annual modulation signa- 
ture exploited by DAMA/Nal; ii) they use different methodological approaches, 
which do not allow any model independent comparison and they have different 
sensitivities to WIMPs; in particular, the number of counts they could expect 
on the basis of the model independent DAMA/Nal result varies from few to 
zero events depending on the models, on the assumptions and on the theoret- 
ical/experimental parameters' values adopted in the calculations; iii) they do 
not make neither correct nor complete comparisons with the DAMA/Nal exper- 
imental result; iv) they use extremely poor statistics; v) they reduce their huge 
measured counting rate of orders of magnitude by various rejection procedures 
claiming for very optimistic rejection powers; vi) their energy scale determina- 
tion and/or energy threshold appear questionable (in the first two cases because 
of the quenching factors values and in the second because of the poor number of 
photoelectrons/keV); etc. 

Finally, in addition, these experiments 8 intrinsically could never reliably claim for 
a signal because of the used approach, as mentioned in §3.2.1. 

5.1.2 ... from indirect searches 

It has been suggested that Dark Matter particles could loose their velocity down to 
a value lower than the escape velocity of a celestial body (Earth, Sun) scattering off 
nuclei and, therefore, remaining trapped in its gravitational field. Subsequently, via 
their annihilation in the celestial bodies or in the Galactic halo they could give rise 
to high energy neutrinos, positrons, antiprotons and gamma's. In principle, the Sun 
could capture WIMPs more effectively than the Earth because of the higher escape 
velocity, but the smaller distance detector - center of the Earth and the "resonant" 
scattering on the heavy nuclei in the Earth (mostly on iron) could compensate this 
effect. 

A possible signature of WIMP annihilation in celestial bodies is given by the pro- 
duced whose interactions in the rock below a detector would give rise to "upgoing" 

"Recent updates of some results, appeared during the publication of this paper, leave the arguments 
discussed in this section unchanched in the essence. 
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muons in the detector itself. The expected \x flux depends on the WIMP annihilation 
rate in the celestial body and on the neutrino energy spectrum produced in the an- 
nihilations. However, several sources of uncertainties are present in similar estimates 
(and, therefore, in the obtained results) such as for example the assumption that a 
"steady state" has been reached in the considered celestial body and the significant 
uncertainty which arises from the estimate and subtraction of the existing competing 
process offered by the atmospheric neutrinos. 

Anyhow, when a model and the related parameters' values are assumed, it is pos- 
sible to estimate the differential flux expected for the secondary neutrinos. According 
to ref. [84] (where the neutralino in the MSSM model has been considered), this flux 
can be written as 

dNu_ ?a ST p( f ) dN f» 

dE u ~ 4nd 2 4- ^ dE v ' [b) 

where Ta is the annihilation rate, d is the distance between the detector and the 
source (e.g. the Earth center or the Sun center), F is the final state of the annihilation 
process, -B^/ are * ne branching ratios of the heavy quarks decays; the dNf v /dE v term 
represents the differential distribution of neutrinos produced by r and by the quarks 
and gluons hadronization and of the subsequent semileptonic decays of the produced 
hadrons. 

Considering, in particular, the and F M , an estimate of the produced neutrino flux 
can be obtained by measuring the up-going muons given by the and T7 M interactions 
with the rock below the detector. Their energy distribution can be written as: 

dN »-N I'™ JF f°° HX f E " dF 1 P (F F ,. x M E ^) dN v 
-dE-- NA ] Ef dEv j dX J E)i dE v P surv (E„E^X) — — , (7) 

where X is the muon range in the rock, da(E^, E'^/dE'^ is the charge current 
cross section for muon production of energy E' from a neutrino of energy E v and 
P S urv{E lll E'^; X) is the survival probability of a muon with E'^ initial energy and E^ 
final energy after crossing a thickness X of rock; Efj 1 is, finally, the energy threshold 
of the detector. The function P surv (E^, E' ; X) obviously account for the muon energy 
loss in the rock. 

As mentioned, the up-going muons produced by atmospheric neutrinos are side 
reactions for the process searched for, however - in principle - they are expected 
to have a flat angular distribution while those induced by WIMPs have a preferred 
impinging direction (e.g. the Sun-laboratory direction or Earth center-laboratory 
direction). 

Model dependent analyses with a similar approach have been carried out by large 
experiments deep underground such as e.g. MACRO and Superkamiokande. It is 
worth to remark that no quantitative comparison can be directly performed between 
the results obtained in direct and indirect searches because it strongly depends on the 
assumptions and on the considered model frameworks. In particular, a comparison 
would always require the calculation and the consideration of all the possible WIMP 
configurations in the given particle model (e.g. for in the allowed parameters space), 
since it does not exist a biunivocal correspondence between the observables in the two 
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kinds of experiments: WIMP-nucleus elastic scattering cross section (direct detection 
case) and flux of muons from neutrinos (indirect detection case). In fact, the counting 
rate in direct search depends on the spin-dependent (SD) and on the spin-independent 
(SI) cross sections of the elastic processes, while the muon flux is connected not only 
to them, but also to the WIMP annihilation cross section. In principle, the three 
cross sections can be correlated, but only when a specific model is adopted and by non 
directly proportional relations. As an example, we report in Fig. 8 the scatter plot 
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Figure 8: Scatter plot for the up-going muon flux from the center of the Earth vs 
neutralino mass. The configurations (MSSM) - here the model dependent constraints 
from the LEP data of 2000 have been used - have been selected by the DAMA annual 
modulation region obtained for the model framework considered in ref . [63] . For details 
see [79]. The solid line is the model dependent upper bound derived from MACRO 
experiment; the one from Superkamiokande is only marginally more stringent. 

for the up-going muon flux from the center of the Earth for a standard Maxwellian 
distribution versus \ mass in MSSM [79]; here the configurations have been selected 
by the DAMA annual modulation region for the particular purely SI model framework 
considered in ref. [63]. As it is evident, the up-going muon flux spans several orders 
of magnitude although the cross section of the DAMA region, allowed in the model 
framework considered there, spans almost one. The solid line in this figure is the 
model dependent upper bound derived from MACRO experiment [85]; the one from 
Superkamiokande is only marginally more stringent, thus it is still compatible with 
the DAMA result even when - as in the quoted ref. [63] - the uncertainties on several 
assumptions and parameters are not yet included (see e.g. [11]). 

As we mentioned at the beginning, the annihilation of the Dark Matter particles 
in the halo could also produce antimatter particles and gamma's. The antimatter 
searches have to be carried out outside the atmosphere, i.e. on balloons or satellites. 
In particular, the WIMP annihilation would result in an excess of antiprotons or of 
positrons up to the WIMP mass with respect to the background arising from other 
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possible sources. Again the estimate and subtraction of such a background together 
with the influence of the Earth and of the galactic magnetic field on these particles 
plays a crucial role on the possibility of a reliable extraction of a signal. However, 
at present some interesting results have been reported in the analysis of the HEAT 
balloon-borne experiment and in some others, as reported in Fig. 9a) and b). In fact, 
an excess of positrons with energy ~ 5 — 20 GeV has been found; it - interpreted in 
terms of WIMP annihilation [87] - gives a result not in conflict with the effect observed 
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Figure 9: Experimental results and theoretical predictions for the indirect searches: 
a) Positron fraction measured by several experiments (Caprice98, AMS, Caprice94, 
HEAT, Clem et al., TS93, MASS89, Golden et al., Muller and Tang, Fanselow et al.; 
see ref. [86] for details). The two lines are two models of background from secondary 
production; an excess of positrons over the background modelling is present in both 
cases, b) Positron fraction measured by HEAT in the interpretation of ref. [87], where 
the excess of positrons with respect to the background modelling is explained in terms 
of neutralino annihilation in the galactic halo, c) Gamma ray energy spectrum of 
the inner Galaxy as measured by EGRET compared with the background modelling 
(lower line). The large excess of gamma's is explained in ref. [88] by a neutralino 
annihilation in the galactic halo (the upper curve is the total contribution). These 
preliminary indications are not in conflict with the DAMA/Nal model independent 
results previously published. 
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by DAMA/Nal. Further results can be expected in future by experiments operating 
in space [89]. 

As regards the possibility to detect 7's from WIMP annihilation in the galactic 
halo, experiments in space are planned. However, at present it is difficult to estimate 
their possibilities considering e.g. the background level, the uncertainties in its reliable 
estimate and subtraction as well as the smallness of the expected signal (even more, if 
a subdominant component would be present) when properly calculated with rescaling 
procedure. However, we mention the analysis of ref. [90] which already suggests the 
presence of a 7 excess from the center of the Galaxy in the EGRET data [91] as 
reported in Fig. 9c). This excess match with a possible WIMP annihilation in the 
galactic halo [90] and is not in conflict with the DAMA/Nal model independent result 
previously published. Other activities are in preparation and will further clarify the 
situation [92]. 

We stress, however, that the specific parameters of a WIMP candidate (mass and 
cross sections), which can be derived from the indirect searches, critically depend on 
several assumptions used in the calculations such as the estimation of the background, 
the halo model, the amount of WIMP in the galactic dark halo, the annihilation 
channels, the transport of charged particle to the Earth, etc.; thus, they have the 
same relative meaning as those obtained in the quest for a candidate in direct search 
approach as described later. 

5.2 Conclusions 

In conclusion, no model independent comparison with the DAMA/Nal effect is avail- 
able. Only few model dependent approaches have been used in the direct search 
approach to claim for a particular model dependent comparison, which appears in 
addition - as discussed above - neither based on solid procedures nor fully correct 
nor complete. On the other hand, the indirect search approaches, which also can offer 
only model dependent comparisons, are either not in contradiction or in substantial 
agreement with the DAMA/Nal observed effect. 

Thus, the interest in the further available DAMA/Nal data is increased; the model 
independent result on the WIMP annual modulation signature from the data of seven 
annual cycles as well as some (of the many possible) model dependent quests for a 
candidate are discussed in the following, reviewing as well the general aspects related 
to the WIMP direct detection. 

6 The DAMA/Nal model independent result on 
the WIMP annual modulation signature from the 
data of seven annual cycles 

As mentioned, the main goal of the DAMA/Nal experiment is the investigation on a 
WIMP component in the galactic halo by exploiting the WIMP annual modulation 
signature introduced in §3.3. The experiment has collected data during seven annual 
cycles. 

The data of each annual cycle have been taken in the same experimental conditions; 
in particular, the Copper box housing the detectors has always been closed and sealed 
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and the detectors have always been in contact only with an atmosphere of HP Nitrogen 
[41, 64]. Moreover, the data taking of each annual cycle has been started before the 
expected minimum of the signal rate (which is roughly around ~ 2 nd December) and 
concluded after the expected maximum (which is roughly around ~ 2 nd June). 

As mentioned, several operational parameters have been regularly acquired with 
the production data, such as the operating temperature, the HP Nitrogen flux into 
the inner Cu box housing the detectors, the pressure of the HP Nitrogen atmosphere 
in the inner Cu box, the environmental Radon from which however the detectors are 
excluded (see above and later) and the hardware rate (including the noise) above 
single photoelectron threshold. Computer controlled processes immediately inform 
the operator during production runs in case one of the parameters goes outside the 
stringent allowed interval of stability values. In addition, the recorded parameters 
values allow a deep analysis and control of possible systematics as performed e.g. in 
ref. [64] and discussed in the following. 

6.1 The evidence 

A model independent investigation of the annual modulation signature has been re- 
alized by exploiting the time behaviour of the residual rates of the single hit events 
in the lowest energy regions over the seven annual cycles (total exposure: 107731 kg 
• day), as previously performed in refs. [63, 64]. These residual rates are calculated 
from the measured event rate after subtracting the constant part (the weighted mean 
of the residuals must obviously be zero over each period): < r^k ~ flatjk >jk- There 
Tijk is the rate in the considered z-th time interval for the j-th detector in the fc-th 
considered energy bin, while flatjk is the rate of the j-th detector in the fc-th energy 
bin averaged over the cycles. The average is made on all the detectors (j index) and 
on all the energy bins in the considered energy interval. 

This model independent approach on the data of the seven annual cycles offers an 
immediate evidence of the presence of an annual modulation of the rate of the single 
hit events in the lowest energy region as shown in Fig. 10, where the time behaviours 
of the (2-4), (2-5) and (2-6) keV single hit residual rates are depicted. They refer to 
4549, 14962, 22455, 16020, 15911, 16608, 17226 kg • day exposures, respectively for 
the DAMA/NaI-1 to -7 running periods 9 . 

In fact, the data favour the presence of a modulated cosine- like behaviour (A- 
cosui(t — to)) at 6.3 a C.L. 10 and their fit for the (2-6) keV larger statistics energy 
interval offers modulation amplitude equal to (0.0200 ± 0.0032) cpd/kg/keV, to = 
(140 ± 22) days and T = ^ = (1.00 ± 0.01) year, all parameters kept free in the fit. 

9 In particular, the DAMA/NaI-5 data have been collected from August 1999 to end of July 2000 
(statistics of 15911 kg • day); then, the DAQ and the electronics have been fully substituted (see 
§4.1). Afterwards, the DAMA/NaI-6 data have been collected from November 2000 to end of July 
2001 (statistics of 16608 kg ■ day), while the DAMA/NaI-7 data have been collected from August 2001 
to July 2002 (statistics of 17226 kg • day), when the data taking with this set-up has been concluded. 

10 It is worth to note that the confidence level given in ref. [63] was instead referred to the particular 
model framework considered there in the quest for a candidate. Here the confidence level refers to 
the model independent effect itself and is calculated on the basis of the residual rate in the (2-6) keV 
energy interval. Applying the same procedure to the residuals given in ref. [63], one gets 4.6 a C.L. 
which is in agreement with the presently quoted value once scaling it by the square root of the ratio 
of the relative exposures. 



28 



The fitting function has been derived from eq. (3) integrated over each time bin. The 
period and phase agree with those expected in the case of a WIMP induced effect (T 
= 1 year and to roughly at ~ 152.5-th day of the year). The \ 2 test on the (2-6) keV 
residual rate in Fig. 10 disfavours the hypothesis of unmodulated behaviour giving a 
probability of 7 • 1CT 4 ( X 2 /d.o.f. = 71/37). We note that, for simplicity, in Fig. 10 
the same time binning already considered in rcf. [63, 64] has been used. The result of 
this approach is similar by choosing other time binnings; moreover, the results given 
in the following are not dependent on time binning at all. 

The residuals given in Fig. 10 have also been fitted, according to the previous 
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Figure 10: Model independent residual rate for single hit events, in the (2-4), (2-5) 
and (2-6) keV energy intervals as a function of the time elapsed since January 1-st of 
the first year of data taking. The experimental points present the errors as vertical 
bars and the associated time bin width as horizontal bars. The superimposed curves 
represent the cosinusoidal functions behaviours expected for a WIMP signal with a 
period equal to 1 year and phase at 2 nd June; the modulation amplitudes have been 
obtained by best fit. See text. The total exposure is 107731 kg • day. 
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procedure, fixing the period at 1 year and the phase at 2 nd June; the best fitted 
modulation amplitudes are: (0.0233 ± 0.0047) cpd/kg/keV for the (2-4) keV energy 
interval, (0.0210 ± 0.0038) cpd/kg/keV for the (2-5) keV energy interval, (0.0192 ± 
0.0031) cpd/kg/keV for the (2-6) keV energy interval, respectively. 

The same data have also been investigated by a Fourier analysis (performed ac- 
cording to ref. [93] including also the treatment of the experimental errors and of the 
time binning), obtaining the result shown in Fig. 11, where a clear peak for a period 
of 1 year is evident. 
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Figure 11: Power spectrum of the measured (2-6) keV single hit residuals calculated 
according to ref. [93], including also the treatment of the experimental errors and of 
the time binning. As it can be seen, the principal mode corresponds to a frequency of 
2.737- 10~ 3 d _1 , that is to a period of ~ 1 year. 



In Fig. 12 the single hit residual rate in a single annual cycle from the total exposure 
of 107731 kg • day is presented for two different energy intervals; as it can be seen the 
modulation is clearly present in the (2-6) keV energy region, while it is absent just 
above. 

Finally, Fig. 13 shows the distributions of the variable s ™~< s ™> , where S m are 
the modulation amplitudes evaluated by a maximum likelihood method [60] for each 
detector, each annual cycle and each considered energy bin (taken there as an example 
equal to 0.25 keV) and a are their errors. The < S m > represent the mean values of the 
modulation amplitudes over the detectors and the annual cycles for each energy bin. 
The left panel of Fig. 13 shows the distribution referred to the region of interest for the 
observed modulation: 2-6 keV, while the right panel includes also the energy region 
just above: 2-14 keV. These distributions allow one to conclude that the individual 
S m values follow a normal distribution, since the variable s ™~< s ™> is distributed as 
a gaussian with an unitary standard deviation. In particular, this demonstrates that 
the modulation amplitudes are statistically well distributed in all the crystals, in all 
the data taking periods and considered energy bins. 
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Figure 12: Single hit residual rate in a single annual cycle from the total exposure 
of 107731 kg • day. The experimental points present the errors as vertical bars and 
the associated time bin width as horizontal bars. The initial time is taken at August 
7 th . Fitting the data with a cosinusoidal function with period of 1 year and phase at 
152.5 days, the following amplitudes are obtained: (0.0195 ± 0.0031) cpd/kg/keV and 
— (0.0009 ± 0.0019) cpd/kg/keV, respectively. Thus, a clear modulation is present in 
the lowest energy region, while it is absent just above. 
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Figure 13: Distributions of the variable ^ (where a is the error associated to 
the S m ) evaluated for each detector, for each annual cycle and each considered energy 
bin: i) in the region of interest for the observed modulation, 2-6 keV (left panel); ii) 
including also the energy region just above, 2-14 keV (right panel). See text. 



6.2 The investigation of possible systematic effects and side 
reactions 



As previously mentioned - to mimic the annual modulation signature a systematic 
effect or side reaction should not only be quantitatively significant, but also able to 
satisfy the six requirements as for a WIMP induced effect (see §3.3); no effect able to 
mimic the signature has been found. A careful investigation of all the known possible 
sources of systematics and side reactions has been regularly carried out by DAMA/Nal 
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Table 3: Modulation amplitudes obtained by fitting the time behaviours of the main 
running parameters including a WIMP-like cosine modulation. These running param- 
eters, acquired with the production data, arc: i) the operating temperature of the 
detectors; ii) the HP Nitrogen flux in the inner Cu box housing the detectors; hi) the 
pressure of the HP Nitrogen atmosphere of the inner Cu box housing the detectors; 
iv) the environmental Radon in the inner part of the barrack from which the detec- 
tors are however excluded (see text); v) the hardware rate above single photoelectron 
threshold. See also the discussion in the whole section. 





DAMA/NaI-5 


DAMA/NaI-6 


DAMA/NaI-7 


Temperature 


-(0.033 ±0.050)°C 


(0.021 ±0.055)°C 


-(0.030 ±0.056)°C 


Flux 


(0.03 ± 0.08) 1/h 


(0.05 ±0.14) 1/h 


(0.07 ±0.14) 1/h 


Pressure 


-(0.6±1.7)10~ 3 mbar 


(0.5±2.5)10~ 3 mbar 


(0.2±2.8)10~ 3 mbar 


Radon 


-(0.09 ±0.17) Bq/m 3 


(0.06 ±0.14) Bq/m 3 


-(0.02 ±0.03) Bq/m 3 


Hardware rate 


(0.10±0.17)10- 2 Hz 


-(0.09 ± 0.19)10- 2 Hz 


-(0.22 ± 0.19)10~ 2 Hz 



and presented at time of each data release [60, 61, 63, 64]. In particular, in ref. [64] 
a detailed discussion has been carried out considering in the quantitative evaluations 
the data of DAMA/NaI-3 and -4. The same analysis is presented in the following 
considering in the quantitative evaluations the data of the DAMA/NaI-5 to -7 running 
periods; it offers a general description of this approach as well. 

First of all the time behaviours of the main running parameters acquired with the 
production data have been investigated. In particular, the modulation amplitudes 
given in Table 3 (see also next subsections) have been obtained for each annual cycle 
when fitting the behaviours including a WIMP-like cosine modulation. As it can be 
seen, all the measured amplitudes are compatible with zero. 

6.2.1 The Radon 

As already discussed elsewhere [41, 64], the detectors are excluded from the environ- 
mental air which contains traces of radioactive Radon gas ( 222 Rn - T ± / 2 = 3.82 days 
- and of 220 Rn - Ti/ 2 = 55 s - isotopes, which belong to the 238 U and 232 Th chains, 
respectively), whose daughters attach themselves to surfaces by various processes. In 
fact: i) the walls, the floor and the top of the installation are insulated by Supronyl 
(permeability: 2 • 10~ n cm 2 /s [94]) and a large flux of HP Nitrogen is released in 
the closed space of that inner barrack; an Oxygen level alarm informs the operator 
before entering it, when necessary; ii) the whole passive shield is sealed in a plexiglas 
box and maintained continuously in HP Nitrogen atmosphere in slight overpressure 
with respect to the environment as well as the upper glove box for calibrating the 
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detectors; iii) the detectors are housed in an inner sealed Cu box also maintained 
continuously in HP Nitrogen atmosphere in slight overpressure with respect to the 
environment; the Cu box can enter in contact only with the upper glove box - during 
calibrations - which is also continuously maintained in HP Nitrogen atmosphere in 
slightly overpressure with respect to the external environment. 

Notwithstanding the above considerations, the Radon in the installation outside 
the plexiglas box, containing the passive shield, is continuously monitored; it is at level 
of sensitivity of the used Radonmeter (see Fig. 14). No modulation of external Radon, 
from which the detectors are anyhow excluded, is observed as quantitatively reported 
in Table 3. 
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Figure 14: Time behaviours of the environmental Radon in the inner part of the bar- 
rack (from which the detectors are however excluded; see text) during the DAMA/Nal- 
5 to -7 running periods, respectively. The measured values are at the level of sensitivity 
of the used radonmeter. 



In Fig. 15 the distribution of the relative variations of the HP Nitrogen flux in 
the inner Cu box housing the detectors and of its pressure as measured during the 
DAMA/NaI-5 to -7 running periods are shown (the typical flux mean value for each 
annual cycle is of order of ~ 260 1/h and the typical overpressure mean value is of 
order of 2 mbar). 

As already reported in rcf. [64], also possible Radon trace in the HP Nitrogen 
atmosphere inside the Cu box has been estimated by searching for the double coinci- 
dences of the gamma-rays (609 and 1120 keV) from 214 Bi Radon daughter, obtaining 
an upper limit on the possible Radon concentration in the Cu box HP Nitrogen atmo- 
sphere: < 4.5 • 10~ 2 Bq/m 3 (90% C.L.); thus, roughly < 4 TO" 4 cpd/kg/keV can be 
expected from this source in the lowest energy bins of interest from the obtained result 
on the WIMP annual modulation signature [64] . This has allowed us to show that even 
an hypothetical, e.g. 10%, modulation of possible Radon in the HP Nitrogen Cu box 
atmosphere would correspond to < 0.2% of the observed modulation amplitude. 
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Figure 15: Distribution of the relative variations of the HP Nitrogen flux in the inner 
Cu box housing the detectors and of its pressure as measured during the DAMA/NaI-5 
to -7 running periods. 



Finally, it is worth to note that, while the possible presence of a sizeable quantity 
of Radon nearby a detector would forbid the investigation of the WIMP annual mod- 
ulation signature (since every Radon variation would induce variation in the measured 
background and the continuous pollution of the exposed surfaces by the non-volatile 
daughters), it cannot mimic the WIMP annual modulation signature in experiments 
such as DAMA/Nal which record the whole energy distribution. In fact, possible pres- 
ence of Radon variation can be easily identified in this case, since it would induce rate 
variation also in other energy regions than the one of interest for the WIMP search, 
that is some of the six requirements of the WIMP annual modulation signature would 
fail. 

In conclusion, no significant effect is possible from the Radon. 
6.2.2 The temperature 

To avoid any significant temperature variation and, in particular, to maintain suitably 
stable the temperature of the electronic devices the installation, where the ~ 100 kg 
Nal(Tl) set-up is operating, is air-conditioned. Moreover, the operating temperature 
of the detectors in the Cu box (stored with the production data) is read out by a probe 
located inside the multi-tons passive shield, whose huge heat capacity assures further 
a relevant stability of the detectors' operating temperature [41, 60, 61, 64, 95]. 

Information for the new DAMA/NaI-5 to 7 running periods can be derived from 
Fig. 16; moreover, no evidence of any operating temperature modulation has been 
observed as quantitatively reported in Table 3. Notwithstanding, to properly evaluate 
the real effect of possible variations of the detectors' operating temperature on the light 
output, we consider - according to the procedure given in ref. [64] - the distribution of 
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Figure 16: Distribution of the relative variations of the operating temperature mea- 
sured during the DAMA/NaI-5 to -7 running periods. 



the root mean square temperature variation within periods with the same calibration 
factors (typically ~ 7 days); this is given in Fig. 17 cumulatively for the three data 
sets. 
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Figure 17: Distribution of the root mean square (r.m.s.) detectors' operating temper- 
ature variation within periods with the same calibration factors (typically ~ 7 days) 
during the DAMA/NaI-5 to -7 running periods. The mean value is 0.05 °C. 

Considering the obtained mean value of the root mean square detectors' operating 
temperature variation: ~ 0.05 °C, and the known value of the slope of the light output 
around its value: < -0.2%/°C, the relative light output variation is < 10~ 4 , which 
corresponds to < 0.5% of the modulation amplitude observed in the lowest energy 
region of the production data, S!^ s , since the counting rate is ~ 1.0 cpd/kg/keV in 
the region of interest [64, 59] and S^ s is ~ 0.02 cpd/kg/keV (see previous section). 
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Thus, any significant effect from the detectors' operating temperature is further 
excluded. As in ref [64], for the sake of completeness, we comment that sizeable 
temperature variations could also induce variations in the electronic noise, in the Radon 
release from the rocks and, therefore, in the environmental background; these specific 
topics will be further analysed in the following, where cautious upper limits on their 
possible effect are given. 

Finally, it is worth to remark that any hypothetical effect induced by temperature 
variations would fail at least some of the six requirements needed to mimic the annual 
modulation signature (such as e.g. the 4 th and the 5 th ). 

In conclusion, all the arguments given above exclude any role of possible effects on 
the observed rate modulation correlated with temperature. 



6.2.3 The noise 

Despite the stringent used noise rejection procedure (see refs. [41, 61, 63, 64] and the 
brief mention in §4.1), the role of possible noise tail in the data after the noise rejection 
has been quantitatively investigated [64] . 

In particular, the hardware rate of each one of the nine detectors above a single pho- 
toelectron, Rnj (j identifies the detector), can be considered; in fact, it is significantly 
determined by the noise. For this purpose the variable Rh = ^j(RHj— < Rnj >) can 
be built (where in our case < Rnj >< 0.25 Hz [41]); its time behaviour during the 
DAMA/NaI-5 to -7 running periods is shown in Fig. 18. 




w 1000 1200 1400 1600 1800 2000 

time (d) 

Figure 18: Time behaviour of the variable Rh = T,j(Rnj — < Raj >), where Rnj is 
the hardware rate of each one of the nine detectors above single photoelectron threshold 
(that is including the noise), j identifies the detector and < Rnj > is the mean value 
of Ruj in the corresponding running period. 

As it can be seen in Fig. 19, the cumulative distribution of Rh for the DAMA/Nal- 
5 to -7 running periods shows a gaussian behaviour with a = 0.5%, value well in 
agreement with that expected on the basis of simple statistical arguments. 

Moreover, by fitting the time behaviour of Rh in the three data taking periods 
- including a WIMP-like modulated term - a modulation amplitude compatible with 
zero: — (0.06±0.11) • 10 -2 Hz, is obtained. From this value the upper limit at 90% C.L. 
on the modulation amplitude can be derived: < 1.3 • 10~ 3 Hz. Since the typical noise 
contribution to the hardware rate of each one of the 9 detectors is ~ 0.10 Hz, the upper 
limit on the noise relative modulation amplitude is given by: ^xQ^jf/ — 1-4 • 10~ 3 
(90% C.L.). Therefore, even in the worst hypothetical case of a 10% contamination 
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Figure 19: Distributions of Rh during the DAMA/NaI-5 to -7 running periods; see 
text. 

of the residual noise - after rejection - in the counting rate, the noise contribution to 
the modulation amplitude in the lowest energy bins would be < 1.4 • 1CP 4 of the total 
counting rate. This means that an hypothetical noise modulation could account at 
maximum for absolute amplitudes of the order of few 10~ 4 cpd/kg/keV, that is <1% 
of the observed annual modulation amplitude [63] . 

In conclusion, there is no evidence for any role of an hypothetical tail of residual 
noise after rejection. 



6.2.4 The efficiencies 

The behaviour of the used efficiencies during the whole data taking periods has even 
been investigated. Their possible time variation depends essentially on the stability 
of the efficiencies related to the previously mentioned acceptance windows, which are 
regularly measured by dedicated calibrations (see e.g. rcf. [41, 61, 63]). 

In particular, we show in Fig. 20 the percentage variations of the efficiency values in 
the (2-8) keV energy interval considering 2 keV bins. They show a gaussian distribution 
with a = 0.5% for DAMA/NaI-5 to -7, cumulatively. Moreover, we have verified that 
the time behaviour of these percentage variations does not show any modulation with 
period and phase expected for a possible WIMP signal. In Table 4 the modulation 
amplitudes of the efficiencies in each energy bin between 2 and 10 keV are reported, 
showing that they are all consistent with zero. In particular, modulation amplitudes - 
considering the three periods together - equal to (0.7± 1.0) • 10~ 3 and (0.1 ± 0.8) • 10~ 3 
are found in the (2-4) keV and (4-6) keV energy bins, respectively; both consistent 
with zero. 

Thus, also the unlikely idea of a possible role played by the efficiency values in the 
observed effect is ruled out. 
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Table 4: Modulation amplitudes obtained by fitting the time behaviour of the effi- 
ciencies including a WIMP-like cosine modulation for the DAMA/NaI-5 to -7 running 
periods. 



Energy 


A 

DAMA/NaI-5 


mplitudc (xlO ' 
DAMA/NaI-6 


! ) 

DAMA/NaI-7 


2-4 kcV 
4-6 kcV 
6-8 kcV 
8-10 keV 


(1.0 ±3.3) 
(1.6 ±2.3) 
(1.0 ± 1.8) 
(0.7 ± 1.3) 


(1.8 ± 1.5) 
(0.7 ± 1.3) 
-(0.1 ±1.0) 
(0.3 ±0.8) 


-(0.4 ± 1.5) 
-(0.9 ± 1.2) 
(0.3 ± 1.0) 
(1.5 ± 1.0) 



6.2.5 The calibration factor 



In long term running conditions, the knowledge of the energy scale is assured by 
periodical calibration with 241 Am source and by continuously monitoring within the 
same production data (grouping them each ~ 7 days) the position and resolution of 
the 210 Pb peak (46.5 keV), mentioned in §4.1 [41, 60, 61, 63, 64]. Although it is highly 
unlikely that a variation of the calibration factor (proportionality factor between the 
area of the recorded pulse and the energy), tdcal, could play any role, according to 
e.g. ref. [64] a quantitative investigation on that point has been carried out. 

For this purpose, the distribution of the relative variations of tdcal - without 
applying any correction - estimated from the position of the 210 Pb peak for all the 
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Figure 20: Distribution of the percentage variations of the efficiency values with the 
respect to their mean values during the DAMA/NaI-5 to -7 running periods; see text. 
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Figure 21: Distribution of the percentage variations of the energy scale factors (tdcal) 
in the DAMA/NaI-5 to -7 running periods without applying any correction; see text. 
The standard deviation is 0.5%. 

9 detectors during the DAMA/NaI-5 to -7 running periods is given in Fig. 21. This 
distribution shows a gaussian behaviour with a < 0.5%. Since the results of the 
routine calibrations are obviously properly taken into account in the data analysis, 
such a result allows to conclude that the energy calibration factor for each detector is 
known with an uncertainty -C 1% within every 7 days interval. 

As discussed also in ref. [64] , the variation of the calibration factor for each detector, 
within each interval of ~ 7 days, would give rise to an additional energy spread (a ca i) 
besides the detector energy resolution (a res ). The total energy spread can be, therefore, 
written as: a = \J a 2 es + <r 2 al ~ a res ■ [1 + | • (f^) 2 ]; clearly the contribution due to 

the calibration factor variation is negligible since | • (f ££i y§) 2 < 7.5 • 10~ 4 20 ^ eV (where 

the adimensional ratio 2 okeV accoun ts for the energy dependence of this limit value). 
This order of magnitude is confirmed by a MonteCarlo calculation, which credits - 
as already reported in ref. [64] a maximum value of the effect of similar variations of 
tdcal on the modulation amplitude equal to 1.6 • 10~ 4 , giving an upper limit < 1% of 
the modulation amplitude measured at very low energy. 

Thus, also the unlikely idea that the calibration factor could play a role can be 
safely ruled out. 

6.2.6 The background 

Similarly as done for the previous data sets (see e.g. ref. [64]), in order to verify the 
absence of any significant background modulation, the energy distribution measured 
during the data taking periods in energy regions not of interest for the WIMP-nucleus 
elastic scattering can be investigated in order to verify if the modulation detected in 
the lowest energy region could be ascribed to a background modulation. In fact, the 
background in the lowest energy region can be essentially due to "Compton" electrons, 
X-rays and/or Auger electrons, muon induced events, etc., which are strictly correlated 
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with the events in the higher energy part of the spectrum. Thus, if a modulation 
with time detected in the lowest energy region would be due to a modulation of the 
background (instead of the possible signal) with time, an equal or larger (sometimes 
much larger) modulation in the higher energy regions should be present. For this 
purpose, we have investigated the rate integrated above 90 keV, R90, as a function of 
the time. 
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Figure 22: Distributions of the percentage variations of R90 with respect to the mean 
values for all the detectors in the DAMA/NaI-5 to -7 running periods; see text. 

In Fig. 22 the distribution of the percentage variations of R90 with respect to the 
mean values for all the nine detectors during the DAMA/NaI-5 to -7 running periods 
is given. They show cumulative gaussian behaviours with a ~ 0.9%, well accounted 
by the statistical spread expected from the used sampling time. This result excludes 
any significant background variation. 

Moreover, fitting the time behaviour of R 90 , a WIMP-like modulation amplitude 
compatible with zero is found in each running period: (0.09±0.32) cpd/kg, (0.06±0.33) 
cpd/kg and -(0.03 ±0.32) cpd/kg for DAMA/NaI-5, DAMA/NaI-6 and DAMA/Nal- 
7, respectively. This excludes the presence of a background modulation in the whole 
energy spectrum at a level much lower than the effect found in the lowest energy 
region; in fact, otherwise - considering the R90 mean values - a modulation amplitude 
of order of tens cpd/kg, that is ~ 100 a far away from the measured value, would be 
present. 

A similar analysis performed on the data collected just above the energy region, 
where the modulation is observed, that is in (6 - 10) keV, gives: -(0.0076 ± 0.0065) 
cpd/kg/keV, (0.0012 ± 0.0059) cpd/kg/keV and (0.0035 ± 0.0058) cpd/kg/kcV for 
the three periods, respectively. 

The results of this subsection also demonstrate that the production data satisfy 
the 4 th requirement of the WIMP annual modulation signature. 

Notwithstanding the results given above already account also for the background 



40 



component due to the neutrons, a further additional independent and cautious analysis 
to estimate their possible contribution has been given in rcf. [64]. In particular, the 
effect of the about 1 m concrete (made from the Gran Sasso rock material) which, as 
known, almost fully surrounds (outside the barrack) the DAMA/Nal passive shield 
- acting as a further neutron moderator - has not been cautiously included in these 
estimates, which are recalled in the following. 

As regards the thermal neutrons, the neutron capture reactions 23 Na(n,j) 24 Na 
and 23 Na{n^) 2im Na (cross section to thermal neutrons equal to 0.10 and 0.43 barn, 
respectively [96]) have been investigated. The capture rate is: ~ 0.2 captures/day/kg 
since the thermal neutron flux has been measured to be 1.08 -10~ 6 neutrons • cm -2 - 
s -1 [38] 11 . Assuming cautiously a 10% modulation of the thermal neutrons flux, the 
corresponding modulation amplitude in the lowest energy region has been calculated 
by MonteCarlo program to be < 10~ 5 cpd/kg/keV, that is < 0.05% of the modulation 
amplitude we found in the lowest energy interval of the production data. In addition, 
a similar contribution cannot anyhow mimic the annual modulation signature since it 
would fail some of the six requirements quoted in §3.3 (such as e.g. the A th and the 
5 th ). 

A similar analysis has also been carried out for the fast neutrons case [64]. From 
the fast neutron flux measured at the Gran Sasso underground laboratory, 0.9 • 10~ 7 
neutrons • cm~ 2 - s _1 [39], the differential counting rate above 2 keV has been estimated 
by MonteCarlo to be ~ 10~ 3 cpd/kg/keV. Therefore, assuming - also in this case - 
cautiously a 10% modulation of the fast neutron flux, the corresponding modulation 
amplitude is < 0.5% of the modulation amplitude found in the lowest energy interval. 
Moreover, also in this case some of the six requirements mentioned above would fail. 

Finally, possible side reactions have been also carefully searched for. The only 
process which has been found as an hypothetical possibility is the muon flux modu- 
lation reported by the MACRO experiment [37]. In fact, MACRO has observed that 
the muon flux shows a nearly sinusoidal time behaviour with one year period and 
maximum in the summer with amplitude of ~ 2 %; this muon flux modulation is 
correlated with the temperature of the atmosphere. A simple calculation to estimate 
the modulation amplitude expected from this process in our set-up has been intro- 
duced in rcf. [64] and is recalled in the following. In fact, the muon flux and 
the yield of neutrons produced by muons measured at the underground Gran Sasso 
National Laboratory (Y) are: $ M ~ 20 muons m _2 d _1 [37] and Y ~ (1 - 7) • 10~ 4 
neutrons per muon per g/cm 2 [97], respectively. Thus, the fast neutron rate produced 
by muons is given by: R n = § u • Y ■ M e ff, where M e ff is the effective mass where 
muon interactions can give rise to events detected in the DAMA set-up. Consequently, 
the annual modulation amplitude in the lowest energy region induced in our experi- 
ment by a muon flux modulation as measured by MACRO [37] can be estimated as: 
S$ — R n - g-e- Jae ■ f single ■ 2%/ (M se t-u P ■ A_B), where g is a geometrical factor, e is 
the detection efficiency for elastic scattering interactions, Jae is the acceptance of the 
considered energy window (E > 2 keV), f single is the "single hit" efficiency and 2% is 
the MACRO measured effect. Since M se t-u P — 100 kg and AE ~ 4 keV, assuming 
the very cautious values g ~ e ~ Jae — f single — 0.5 and M e ff = 15 t, one obtains: 

11 Consistent upper limit on the thermal neutron flux have been obtained with the ~ 100 kg DAMA 
Nal(Tl) set-up considering these same capture reactions [41] . 
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Srn < (1 - 7) • 1CT 5 cpd/kg/keV, that is < 0.3% of the modulation amplitude we 
observe [63]. We stress that - in addition - the latter value has been overestimated 
of orders of magnitude both because of the extremely cautious values assumed in the 
calculation and, as mentioned, of the omission of the effect of the ~ 1 m concrete 
neutron moderation. Finally, we remark that not only the modulation of the muon 
flux observed by MACRO would give rise in our set-up to a quantitatively negligible 
effect, but - in addition - some of the six requirements necessary to mimic the annual 
modulation signature (such as e.g. the 4 th and the 5 th ) would fail. Therefore, it can 
be safely ignored. 

Just for the sake of completeness, we remind that the contribution of solar neutri- 
nos, whose flux is also expected to be modulated, is many orders of magnitude lower 
than the measured rate [98]. 

In conclusion, the results presented in this section demonstrate that the production 
data satisfy the A th requirement of a WIMP induced effect and - at the same time 
- exclude that the annual modulation observed in the lowest energy region could be 
ascribed to modulation of any kind of possible background. 

6.3 Conclusions on the DAMA/Nal model independent result 

No modulation has been found in any of the considered possible source of systematics; 
thus, upper limits (90% C.L.) on the possible contributions to a modulated ampli- 
tude have been calculated and are summarized in Table 5. In particular, they cannot 
account for the measured modulation because quantitatively not relevant and, as dis- 
cussed in details already in ref. [64], none of them is able to mimic the observed effect; 
in fact, none can satisfy all the above mentioned peculiarities of the signature. 

The quantitative investigations discussed above offer a complete analysis of known 
sources of possible systematic effects. We can conclude that a relative systematic error, 
affecting the energy spectrum, of order of < 10~ 3 is credited by these investigations, 
while the results on the analysis of Rg exclude the presence of a possible background 
modulation even at more stringent level. Furthermore, no systematic effect or side 
reaction able to mimic a WIMP induced effect, that is to be not only quantitatively 
significant, but also able to satisfy all the many requirements of the signature (see 
§3.3), has been found. 

In conclusion, the presence of an annual modulation in the residual rate in the 
lowest energy interval (2 - 6) keV with all the features expected for a WIMP component 
in the galactic halo is supported by the data at 6.3 a C.L. 

7 Corollary results: quests for a candidate particle 
in some model frameworks 

On the basis of the previous main result, a corollary investigation can also be pur- 
sued on the nature and coupling of a WIMP candidate. This latter investigation is 
instead model dependent and - considering the large uncertainties which exist on the 
astrophysical, nuclear and particle physics assumptions and parameters needed in the 
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Tabic 5: Summary of the results obtained by investigating possible sources of system- 
atics or of side reactions in the data of the DAMA/NaI-5 to 7 running periods. None 
able to give a modulation amplitude different from zero has been found; thus cautious 
upper limits (90% C.L.) on the possible contributions to a modulation amplitude have 
been calculated and are shown here in terms of the measured model independent mod- 
ulation amplitude, S^ s , we have observed (see §6.1). As it can be seen none (nor their 
cumulative) effect could account for the measured modulation; moreover, as discussed 
in details already in rcf.[64], none of them could mimic the signature. 



Source 


Main comment 


Cautious upper limit 
(90%C.L.) 


Radon 


Sealed Cu Box in 
HP Nitrogen atmosphere 


< 0.2%S% S 


Temperature 


Air conditioning 


< 0.5%S% S 


Noise 


Efficient rejection 


< \%S^ S 


Energy scale 


Routine 
+ intrinsic calibrations 


< \%S^ S 


Efficiencies 


Mainly routine measurements 


< i%st 


Background 


No modulation observed 
above 6 keV; this limit 
includes possible effect 
of thermal and fast neutrons 


< 0.5%S% S 


Side reactions 


From muon flux variation 
measured by MACRO 




In addition: no effect can mimic the signature 



calculations - has no general meaning (as it is also the case of exclusion plots and of 
the WIMP parameters evaluated in indirect search experiments). Thus, it should be 
handled in the most general way as we have preliminarily pointed out with time in the 
past [60, 61, 62, 63, 64, 65, 66, 11] and we will again show in the following sections. 
The results we will discuss are, of course, not exhaustive and many other different 
allowed regions can be obtained by varying the assumptions within the many possible 
model frameworks, which at present level of knowledge cannot be disentangled (e.g. 
open questions: i) which is the right nature for the WIMP particle; ii) which is its 
right couplings with ordinary matter; iii) which are the right form factors and related 
parameters for each target nucleus; iv) which is the right spin factor for each target 
nucleus; v) which are the right scaling laws; vi) which is the right halo model and re- 
lated parameters; vii) which are the right values of the experimental parameters which 
we can determine only with associated uncertainties; etc.). The situation is analogous 
for results presented in form of exclusion plots, which also have not an "universal" 
meaning, as well as for the results of WIMP specific parameters inferred by indirect 
searches. 
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7.1 Estimates of WIMP direct detection rates 



In the following the main bases necessary to perform model dependent analyses in 
this field are given. The approximations and intrinsic uncertainties can be inferred 
straight-forward . 



7.1.1 WIMP-nucleus elastic scattering 

The studied process is the WIMP-nucleus elastic scattering and the measured quantity 
is the recoil energy. In the general form the differential energy distribution of the recoil 
nuclei can be calculated [46, 99] by means of the differential cross section of the WIMP- 
nucleus elastic processes 



da / da \ / da \ 

dE^^ ER) -{d^) SI + {dE- R ) : 

2G 2 F m N 



SD 

{ [Zg p + (A- Z)g n } 2 F 2 SI (E R ) + 8A 2 J {J + 1)F 2 D (E R )}, (8) 



where: Gf is the Fermi coupling constant; toat is the nucleus mass; v is the WIMP 
velocity in the laboratory frame; E R = m 2 VN v 2 (l — cos9*)/m,M (with mwN WIMP- 
nucleus reduced mass and 9* scattering angle in the WIMP-nucleus cm. frame) is 
the recoil energy; Z is the nuclear charge and A is the atomic number; g pn are the 
effective WIMP-nucleon coupling strengths for SI interactions; A 2 J (J + 1) is a spin 
factor. Moreover, FgjiEp) is the SI form factor (see later), while Fg D (E R ) is the SD 
form factor (see later) for which an universal formulation is not possible since in this 
case the internal degrees of the WIMP particle model (e.g. supersymmetry in case of 
neutralino) cannot be completely decoupled from the nuclear ones. It is worth to note 
that this adds significant uncertainty in the model dependent results. 

Furthermore, it can be demonstrated that A — a p <s p>+^ n<Sn> w ith J nuclear 
spin, with o Pi „ effective WIMP-nucleon coupling strengths for SD interaction and with 
< Sp. n > mean values of the nucleon spins in the nucleus. Therefore, the differential 
cross section and, consequently, the expected energy distribution depend on the WIMP 
mass and on four unknown parameters of the theory: g p ^ n and a Pi „. 

The total cross section for WIMP-nucleus elastic scattering can be obtained by 

2m 2 v 2 

integrating equation (8) over E R up to E R . max = m ™" v : 



ct(v) = / - nr (v,E R )dE R ^-G 2 F m 2 WN {[Zg p + (A-Z)g n } 2 G S i(v) + 

JO dE R 7T 

+8 ^±1 [a p <S P > +a n < S n >f G SD (v)}. (9) 

Here G S i(v) = ^ R ' max F 2 j(E R )dE R ; G SD (v) can be derived straightforward. 

The standard point-like cross section can be evaluated in the limit v — > (that is 
in the limit Gsi{v) and Gsd{v) — ► 1). Knowing that < 5 Pi „ >= J = 1/2 for single 
nucleon, the SI and SD point-like cross sections on proton and on neutron can be 
written as: 

SI _ ^ /~<2 2 2 SD _ 32 3 2 2 2 /i n \ 

a p,n — ~ { - T F m W(p,n)9p,n a p,n ~ ~ F m W(p,n) a p,ny \ W ) 
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where mwp — mwn are the WIMP-nucleon reduced masses. 

As far as regards the SI case, the first term within squared brackets in eq. (9) can 
be arranged in the form 



[* fc + (il -* to J»_(&±&)"[l-l£i 



A 



A 2 =g 2 -A 2 . (11) 



Considering ^ nearly constant for the nuclei typically used in direct searches for Dark 
Matter particles, the coupling term g is generally assumed - in a first approximation 
- as independent on the used target nucleus. Under this assumption, the nuclear 
parameters can be decoupled from the particle parameters and a generalized SI WIMP- 
nucleon cross section: a si — ■^G 2 F m 2 v , p g 2 , can be conveniently introduced. 

As far as regards the SD couplings, let us now introduce the useful notations [65] 

a = JaJTa 2 , tgO = — , a SD = — ^G^m^a 2 , (12) 

where <tsd is a suitable SD WIMP-nucleon cross section. The SD cross sections on 
proton and neutron can be, then, written as: 

<?p D = <?sd ■ cos 2 6 a s n D = asD ■ sin 2 9. (13) 

In conclusion, equation (8) can be re- written in terms of a si, <Jsd and as: 

' hT (v,E R ) = -^-,.E(E R ), (14) 



dE R 2m w P u 
with 

X(E R ) = {A 2 <j S iF 2 I (E R ) + 

+ Y—f^°SD [< S p > cos#+ < S n > sine} 2 F 2 D (E R )}. (15) 

The mixing angle 9 is defined in the [0, n) interval; in particular, 9 values in the second 
sector account for a p and a n with different signs. As it can be noted from its definition 
[99], Fg D (E R ) depends on a p and a n only through their ratio and, consequently, 
depends on 6, but it does not depend on a. 

Finally, setting the local WIMP density, pw = £po, where po is the local halo 
density and £ 12 (£ < 1) is the fractional amount of local WIMP density, and the 
WIMP mass, mw, one can write the energy distribution of the recoil rate (R) in the 
form 

dR AT p w [ Vma * da 



dE R 



= Nt P^L [ ma * ^L {v ,E R )vf(v)dv = 



Po ■ m N 



^-&{E R )I(E R ), (16) 

2m w ■ m 2 Wp 



12 Pay attention that in ref. [61, 63, 65, 66] the same symbol indicates instead a different quantity: 
£ = Pw /(0.3GeVcm- 3 ). 
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where: Nt is the number of target nuclei and I (En) = J^ maX f E \ dv^p- with f(v) 



WIMP velocity providing E R recoil energy; v max is the maximal WIMP velocity in 
the halo evaluated in the Earth frame. 

The differential distribution of the detected energy, E det , for a multiple- nuclei de- 
tector (as e.g. the Nal(Tl)) can be easily derived: 



where q x is the quenching factor for the x recoiling nucleus and K(Ea et , E') takes into 
account the response and energy resolution of the detector; generally it has a gaussian 
behaviour. 

It is worth to remark, as it can be inferred by eq. (8), that only nuclei with spin 
different from zero are sensitive to WIMPs with both SI and SD couplings. This is the 
case of the 23 Na and 127 I nuclei, odd-nuclei with an unpaired proton, constituents of 
the DAMA/Nal detectors. Thus, the purely SI coupling scenario widely considered in 
this field represents only a particular case of the more general framework of a WIMP 
candidate with both mixed SI and SD couplings. Therefore, in the following analyses, 
we will consider some of the possible scenarios for the mixed SI and SD couplings and, 
then, also the sub-cases of pure SI and pure SD couplings. 

7.1.2 WIMPs with preferred inelastic scattering 

It has been suggested [27] also the possibility that the annual modulation of the low en- 
ergy rate observed by DAMA/Nal could be induced by WIMPs with preferred inelastic 
scattering: relic particles that cannot scatter elastically off nuclei. As discussed in ref. 
[27] , the inelastic Dark Matter could arise from a massive complex scalar split into two 
approximately degenerate real scalars or from a Dirac fermion split into two approx- 
imately degenerate Majorana fermions, namely \+ an d X-, with a 5 mass splitting. 
In particular, a specific model featuring a real component of the sneutrino, in which 
the mass splitting naturally arises, has been given in ref. [27] and mentioned here in 
§2. The detailed discussion of the theoretical arguments on such inelastic Dark Matter 
can be found in ref. [27]. In particular, there has been shown that for the \- inelastic 
scattering on target nuclei a kinematical constraint exists which favours heavy nuclei 
(such as 127 I) with respect to lighter ones (such as e.g. ™ at Ge) as target-detectors me- 
dia. In fact, x- can on ly inelastically scatter by transitioning to x+ (slightly heavier 
state than x_) and this process can occur only if the X- velocity, v, is larger than: 



This kinematical constraint becomes increasingly severe as the nucleus mass, m^, is 
decreased [27]. For example, if 5 > 100 keV, a signal rate measured e.g. in Iodine will 
be a factor about 10 or more higher than that measured in Germanium [27]. Moreover, 
this model scenario implies some peculiar features when exploiting the WIMP annual 



WIMP velocity distribution in the Earth 





(17) 




(18) 
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modulation signature [53]; in fact - with respect to the case of WIMP elastically 
scattering - it would give rise to an enhanced modulated component, S m , with respect 
to the unmodulated one, 5*0, and to largely different behaviours with energy for both 
So and S m (both show a higher mean value) [27]. 

The preferred inelastic Dark Matter scenario [27] offers further possible model 
frameworks and has also the merit to naturally recover the sneutrino as a WIMP 
candidate (see e.g. §7.1.2). 

The differential energy distribution of the recoil nuclei in the case of inelastic pro- 
cesses can be calculated by means of the differential cross section of the WIMP-nuclcus 
inelastic processes: 



= [Zgp + {A _ Z) 9n ?F 2 SI {f) • V 1 - ^ (W) 

where dQ* is the differential solid angle in the WIMP-nucleus cm. frame; q 2 is the 
squared three-momentum transfer. 

In the inelastic process the recoil energy depends on the scatter angle, 0*, in the 
cm. frame according to: 



2 m 2 sitl 2 cos9* 

E R = zm WN v m V (20) 

rriN 2 

Thus, we can write: 



dE R ^<^Jl- V \.^. (21) 

WLM V V Z 47T 

From cq. (19) and (21) we derive the differential cross section as a function of the 
recoil energy, E R , and the WIMP velocity, v. 

(«> En) = [Z 9p + (A- Z)g n f Fg I (E R ). (22) 

Here we apply the relation q 2 = 2mj^E R . 

The minimal WIMP velocity, v m i n (E R ), providing E R recoil energy in the inelastic 
process is: 



m N E R ( m WN S\ 

and it is always > v t hr- 

Finally, one can write the energy distribution of the recoil rate (R) in the form 

dR pw f Vma:c da 



Nt £»L f ma * ^- {v , ER)vf{v)dv = 

m W Jv^(En) dE R 



dE R m w J Vmin (E R ) dE R 

Nt / '^ ■ A 2 ta p F 2 SI (E R ) ■ I(E R ). (24) 
2m w ■ mfap 

Moreover, as derived in the case discussed in the previous subsection, also in 
the present case a generalized SI point-like WIMP-nucleon cross section: a p — 
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^G^m^pg 2 , can be denned. Finally, the extension of formula (24) e.g. to detec- 
tors with multiple nuclei can be easily derived. 

In this scenario the modulated and the unmodulated components of the signal are 
function of £<7 P , mw and S. 



7.1.3 The halo models 

As discussed above, the expected counting rate for the WIMP elastic scattering de- 
pends on the local WIMP density, pw, and on the WIMP velocity distribution, f(v), 
at Earth's position. The experimental observations regarding the dark halo of our 
Galaxy do not allow to get information on them without introducing a model for the 
Galaxy matter density. An extensive discussion about models has been reported in 
ref. [11]. Here we present a brief introduction on this argument both to allow the 
reader to understand the complexity of this aspect and in the light of the results given 
in following subsections on the discussed quests. 

Important information on the dark halo in the Galaxy can be derived from mea- 
surements of the rotational velocities of objects bounded in the gravitational galac- 
tic field. In fact, the following relation between the rotational velocity of an ob- 
ject placed at distance r = |f| from the center of the Galaxy and the total mass, 
M to t(r) = f r , <r d 3 r' p to t(r'), inside the radius r can be obtained from the virial theo- 
rem: 

vUr) _ GM^M (25) 

where G is the Newton's constant. The total mass density, ptot(r), can be expressed 
as the sum of the mass density of the dark halo, Pdm(j), and the mass density of the 
visible component, p v i s (r) that constitutes the bulge and the disk of the Galaxy. The 
gravitational potential ^f(r) is related to p to t through the Poisson's equation: 

V 2 * = ~4nGp tot . (26) 

The halo density profile Pdm(?) can also be expressed in term of the distribution 
function of the WIMP in the six-dimensional phase space F(r, v): 

Pnu(f) = J F(r,v)d 3 v (27) 

where r and v represent the position and velocity vectors in the rest frame of the 
Galaxy respectively. 

Inverting eq. (27) and taking into account observational data it is, in principle, 
possible to calculate F(r, v) and, then, the WIMP velocity distribution function at the 
Earth position in the Galaxy: 

f(v) = F(R ,v) (28) 

where Ro — 8.5 kpc corresponds to the Earth distance from the center of the Galaxy 
along the direction of the galactic plane (Ro = (Ro,0, 0)). The inversion of eq. (27) 
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presents degeneracy problem that can be solved only requiring some degree of symme- 
try for the system. The velocity distribution function represents therefore an impor- 
tant source of uncertainties in the evaluation of expected counting rate for the WIMP 
component in the galactic halo. 

As well as the velocity distribution function, two other quantities are of great 
importance to estimate the expected counting rate: the WIMP local velocity, vq = 
Vrot(Ro)i an d the local halo density /?o = Pdm(Rq) that appears as a multiplicative 
factor in the formula giving the expected counting rate. Since we are interested in 
the evaluation of local and rotational velocity of the WIMP at distance r — Ro, it is 
not necessary for our purpose a detailed description of the inner part of the Galaxy 
(r <C Ro) where disk and bulge are dominant. Thus, it is possible to consider the 
bulge as a spherical density distribution with relevant contribution up to about 1 kpc 
and truncated at r ~ 2 kpc from the center of Galaxy, and the disk as an exponential 
distribution which (according to most of the models) decreases up to 4 kpc where it 
can safely be neglected. Therefore, it is generally assumed that for r > R the dark 
matter is the dominant component. 

The contribution of the visible matter has been considered in the calculation of the 
WIMP local velocity: 

vl = v 2 rot (R ) = [M vis (Ro) + M halo (Ro)] . (29) 

A maximal halo, p nax , occurs when M vis (R ) <§C Mhaio(Ro)] m this case the contri- 
bution to the rotational velocity is due to the halo; on the other hand, when for M, 
the maximum value compatible with observations is considered, a minimal halo, p% 
occurs and only a fraction of v is supported by the dark halo. 

The WIMP halo can be represented as a collisionless gas of particles whose distri- 
bution function satisfy the Boltzmann equation [100]. In general case the Boltzmann 
equation cannot be solved without reducing the complexity of the system. It is pos- 
sible to consider models based on the Jeans or on the virial equations that describe 
a wide range of systems, but one cannot be sure that these models describe systems 
with realizable equilibrium configuration [100]. 

The dark halo model widely used in the calculations carried out in the WIMP direct 
detection approaches is the simple isothermal sphere that corresponds to a spherical 
infinite system with a flat rotational curve. The halo density profile is: 

PDM(r) = £^ (30) 
corresponding to the following potential: 

*o(r) = -| log(r 2 ). (31) 

In this case, when a maximal halo density is considered, the WIMP velocity distribu- 
tion is the Maxwell function: 

f(v) = Acxp (-^-) (32) 



VIS 

min 



V 2v, 



2 

rms 
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where N is the normalization constant. The mean square velocity results: v 2 ms = 
(3/2)t> 2 ; this relation descends from the hypothesis of an halo formed by particles in 
hydrostatic equilibrium with an isotropic velocity distribution. Despite the simplicity 
of this model has favoured its wide use in the calculation of expected rate of WIMP- 
nucleus interaction, it doesn't match with astrophysical observations regarding the 
sphericity of the halo and the absence of rotation, the flatness of the rotational curve 
and the isotropy of the dispersion tensor, and it presents unphysical behavior: the 
density profile, in fact, has a singularity in the origin and implies a total infinite mass 
of the halo unless introducing some cut-off at large radii. 

In the ref. [11] the analysis of the first 4 DAMA/Nal annual cycles in a particu- 
lar case for a SI coupled WIMP candidate has been extended by considering a large 
number of self-consistent galactic halo models, in which the variation of the veloc- 
ity distribution function is originated from the change of the halo density profile or 
of the potential. The different models have been classified in 4 classes according to 
the symmetry properties of the density profile or of the gravitational potential and 
of the velocity distribution function. The same strategy has been followed to obtain 
the new cumulative results given later. The considered halo model classes correspond 
to: spherically symmetric matter density with isotropic velocity dispersion (A); spher- 
ically symmetric matter density with non- isotropic velocity dispersion (B); axisym- 
metric models (C); triaxial models (D). The models are summarized in Table 6 where, 
according to ref. [11], arc identified by a label. 

We will present briefly in the following these models since they will be considered 
in the new results on the quest for possible candidate particle given in the following 
subsections. For a detailed discussion refers to the ref. [11]. 

/. Spherical halo models with isotropic velocity dispersion (A) 

The first class groups models with spherical density profile; for these models p{r) = 
p{r) and f(v) = f(v). The first type of model is a generalization of the spherical 
isothermal sphere in which a core radius R c is introduced. The density profile becomes 
(model Al): 

vl ml + r 2 

PM(r) = 4rfMT^' (33) 

with, in case of maximal halo, corresponding potential: 

*o(r) = -|log( J R c 2 +r 2 ). (34) 

In the limit R c — > the profile (30) and the potential (31) is obtained. These models 
are also named logarithmic because of the analytic form of the potential. 

A second class of spherical models are defined by the following matter density 
profile (models A2 and A3): 

_ p* a RP 3^ + r 2 (l-/3) 
Pdm (r) - {Rl + r2){p+A)/2 , (35) 
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Tabic 6: Summary of the consistent halo models considered in the analysis of ref. [11] 
and in the following. The labels in the first column identify the models. In the third 
column the values of the related considered parameters are reported [11]; Other choices 
are also possible as well as other halo models. In the last column references to the 
corresponding equations in the text are listed. The models of the Class C have also 
been considered including possible co-rotation and counter-rotation of the dark halo 
(see eq. (44).) 



Class A: spherical pr>M, isotropic velocity dispersion 


eq. 


AO 


Isothermal Sphere 




(30) 


Al 


Evans' logarithmic [101] 


R c — 5 kpc 


(33) 


A2 


Evans' power-law [102] 


R c = 16 kpc, = 0.7 


(35) 


A3 


Evans' power-law [102] 


R c = 2 kpc, p = -0.1 


(35) 


A4 


Jaffe [103] 


a = 1, p = 4, 7 = 2, a = 160 kpc 


(37) 


A5 


NFW [104] 


a = 1, P = 3, 7 = 1, a = 20 kpc 


(37) 


A6 


Moore et al. [105] 


a = 1.5, P = 3, 7 = 1.5, a = 28 kpc 


(37) 


A7 


Kravtsov et al. [106] 


a = 2, p = 3, 7 = 0.4, a = 10 kpc 


(37) 


Class B: spherical /)dm, non isotropic velocity dispersion 




(Osipkov Merrit, p = 0.4) 






Bl 


Evans' logarithmic 


R c = 5 kpc 


(33)(39) 


B2 


Evans' power-law 


R c = 16 kpc, P = 0.7 


(35)(39) 


B3 


Evans' power-law 


R c = 2 kpc, P = -0.1 


(35)(39) 


B4 


Jaffe 


a = 1, p = 4, 7 = 2, a = 160 kpc 


(37)(39) 


B5 


NFW 


a = 1, P = 3, 7 = 1, a = 20 kpc 


(37)(39) 


B6 


Moore et al. 


a = 1.5, P = 3, 7 = 1.5, a = 28 kpc 


(37)(39) 


B7 


Kravtsov et al. 


a = 2, p = 3, 7 = 0.4, a = 10 kpc 


(37) (39) 


Class C: Axisymmetric pdm 




CI 


Evans' logarithmic 


Rc = 0, q = l/y/2 


(40)(41) 


C2 


Evans' logarithmic 


R c = 5 kpc, q = 1/V2 


(40)(41) 


C3 


Evans' power-law 


R c = 16 kpc, q = 0.95, P = 0.9 


(42)(43) 


C4 


Evans' power-law 


R c = 2 kpc, q = l/\/2, P = -0.1 


(42)(43) 


Class D: Triaxial p DM [107] (q = 0.8, 


p = 0.9) 




Dl 


Earth on maj. axis, rad. anis. 


S = -1.78 


(45)(46) 


D2 


Earth on maj. axis, tang. anis. 


S = 16 


(45) (46) 


D3 


Earth on interm. axis, rad. anis. 


(5 = -1.78 


(45) (46) 


D4 


Earth on interm. axis, tang. anis. 


S = 16 


(45)(46) 



and potential for a maximal halo: 

*° (r) = {Rl + (/ ^ 0) - (36) 

We will refer to these models as power-law halo models. They represent the spherical 
limit of the more general axisymmetric model discussed later. When the parameter 
P — > 0, the logarithmic models are obtained. 

The last family of spherical models is described by the matter density distribution 
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(models A4 - A7): 

PDM y) = Po 



l + (i?o/a) Ql( ^ 7)/a 



(37) 



1 + (r/a) a 

The different choice of the parameters: a, (3, 7 and a, used in the calculations given 
later are reported in Table 6; other choices are possible. The density profile of these 
models, except for the Jaffe case, has been obtained from numerical simulations of 
Galaxy evolution. 

II. Spherical halo models with non-isotropic velocity dispersion (B) 

These models have been studied in the simple case in which the velocity distribution 
function depends on the two integrals of motion energy and angular momentum vector 
(L = \L\) only through the so called Osipkov-Merrit variable [100, 108]: 

Q = *-^2> (38) 

where the parameter r a appears in the definition of the /3 0) the degree of anisotropy 
of the velocity dispersion tensor on the Earth's position [108]: 

A) = 1-S = "^T- ( 39 ) 
Vr Rq + r 2 a 

In this definition the velocity is expressed in spherical coordinates and v$ = ve 7^ Vr 
(con v? =< vf > - < Vi > 2 ,i = <fi,6,r). 

The considered models are the same as in the isotropic case and the velocity dis- 
tribution function has been calculated introducing the Osipkov-Merrit term in the 
equations. The degree of anisotropy of the models depends on the /?o values; for 
/?o — > 1, or v^ 2 — Tv 2 , the distribution function becomes isotropic. 

777. Axisymmetric models ( C) 

In these models the velocity distribution depends in general at least on the energy e 
and on the component L z of the angular momentum along the axis of symmetry. The 
velocity distribution can be written as the sum of an even and an odd contribution 
with respect to L z . It can be shown [101, 102] that the puu depends only on the 
even part and the velocity distribution can be calculated up to an arbitrary odd part. 
The axisymmetric generalizations of the Evans' logarithmic and power-law models 
have been considered in ref. [11]. For these models the velocity distribution has been 
calculated analytically by Evans [101, 102] and corresponds to a maximal halo. The 
axisymmetric logarithmic potential (models CI and C2) is: 



*oOR,*) = -f log (i? 2 + i? 2 + £), 



(40) 



where R — (x 2 + y 2 ), is the radial coordinate along the galactic plane and R c is the 
core radius; q is the flatness parameter. The corresponding matter density distribution 
results: 

v 2 (2 q 2 + l)R 2 c +R 2 + (2- q - 2 )z 2 
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If an asymptotically non-fiat rotational curve is considered, the axisymmetric 
power-law potential is obtained [102] (models C3 and (74): 

^^wrwfk^w ( ^ 0) - (42) 

with the distribution function: 

_ [M a RP (2g 2 + l)Rl + (1 - /jg 2 )^ + [2 - <T 2 (1 + /3)]z 2 

^M(ii, Z) - {Rl+R 2 +z 2 q -2 ){P+ ,),2 • ( 43 ) 

The related velocity distribution functions for these two cases can be found in ref. 
[11]. 



IV. Co-rotating and counter-rotating halo models 



In the case of axisymmetric models it is possible to include an halo rotation con- 
sidering that the velocity distribution function is known up to an arbitrary odd com- 
ponent. An odd component of velocity distribution function can easily be defined 
starting from an even solution. The velocity distribution function, linear combination 
of even and odd function, is able to describe an halo configuration where a particle 
population moves clockwise around the axis of symmetry and a population moves in 
opposite sense. In this case the velocity distribution can be written as [11]: 

F{e,L z )= V F right (e,L z ) + {l-ri)F left (e,L z ). (44) 

The r\ parameter ranges from 1 (maximal co-rotation) to (maximal counter-rotation) 
and it is related to the dimensionless spin parameter A of the Galaxy by: A = 0.36|?7 — 
0.5| [109]. 

Considering the limit A < 0.05 obtained from numerical work on Galaxy formation 
[110], n can range in the interval 0.36 < r\ < 0.64. For the Evans' axisymmetric models 
of class C we have also considered possible co-rotation and counter-rotation of the 
halo assuming i] = 0.36 and n = 0.64. 



V. Triaxial models (D) 

The models, belonging to this class, arise from the triaxial potential discussed in 
ref. [107]: 

I / y2 z 2\ 

V (x,y,z) = --v 2 log\x 2 + ^ + — J, (45) 
This potential, in the case of a maximal halo, corresponds to the density profile: 



Ax 2 +By 2 /p 2 + Cz 2 /q 2 
4ttG (x 2 +y 2 /p 2 + z 2 /q 2 f 



PDM(x,y,z)= —————— (46) 



where A = (p- 2 +q- 2 -l), B = (l + q- 2 -p- 2 ) and C = (l+p- 2 -q- 2 ). Inref.[107], 
the velocity distribution is approximated by a Gaussian with v 2 , v 2 , vi depending on 
vq, 5, p and q parameters an on the Earth position: i) Earth on the major axis of 
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the equipotential ellipsoid (models Dl and D2); ii) Earth on the intermediate axis 
(Models D3 and DA). In these cases the free parameter 5 appears; this parameter, 
in the limit of spherical halo (p = q = 1), measures the degree of anisotropy of the 

-2 

velocity dispersion tensor: = 



VI. Constraining the models 

The parameters of each halo model, given above, have been chosen in ref. [11] 
taking into account the available observational data. Anyhow, information on galactic 
dark halo can be obtained only in indirect way [111, 112] and considering hypotheses 
on its form and characteristic. 

The allowed range of values for the WIMP local velocity has been estimated there 
considering the information coming from the rotational curve of our Galaxy. The 
considered interval is: 

v = (220 ±50) kms" 1 (90% C.L.), (47) 

that conservatively relies on purely dynamical observations [113]. Similar estimates of 
the vq central value with smaller uncertainty have been obtained studying the proper 
motion of nearby stars in the hypothesis of circular orbit of these objects [114]. In the 
analyses given in the following, for simplicity, we have considered the three represen- 
tative values of local velocity: vo — 170, 220, 270 km/s. 

Table 7: Allowed intervals of po for the halo models quoted in Table 6 as evaluated in 
ref. [11]. The p nax and p™ m values (in GeV cm~ 3 ) are here used in the quests for 
a candidate particle (see later) for the class A and B, while only the case of p™ ax is 
used for the class C and D. 





vq = 170 km s 1 


v Q = 220 km s" 1 


v Q = 270 km s" 1 


Model 


Po 


'max 

Po 


Po 


„max 

Po 


Po 


Po 


AO 


0.18 


0.28 


0.30 


0.47 


0.45 


0.71 


Al , Bl 


0.20 


0.42 


0.34 


0.71 


0.62 


1.07 


A2 , B2 


0.24 


0.53 


0.41 


0.89 


0.97 


1.33 


A3 , B3 


0.17 


0.35 


0.29 


0.59 


0.52 


0.88 


A4 , B4 


0.26 


0.27 


0.44 


0.45 


0.66 


0.67 


A5 , B5 


0.20 


0.44 


0.33 


0.74 


0.66 


1.11 


A6 , B6 


0.22 


0.39 


0.37 


0.65 


0.57 


0.98 


A7 , B7 


0.32 


0.54 


0.54 


0.91 


0.82 


1.37 


CI 


0.36 


0.56 


0.60 


0.94 


0.91 


1.42 


C2 


0.34 


0.67 


0.56 


1.11 


0.98 


1.68 


C3 


0.30 


0.66 


0.50 


1.10 


0.97 


1.66 


C4 


0.32 


0.65 


0.54 


1.09 


0.96 


1.64 


Dl , D2 


0.32 


0.50 


0.54 


0.84 


0.81 


1.27 


D3 , D4 


0.19 


0.30 


0.32 


0.51 


0.49 


0.76 
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Moreover, in the analyses given in the following, for each model - after fixing 
the local velocity - the local density po has been assumed to be equal either to the 
pmin or ^ ne p-max vame as obtained in ref. [11] when imposing the following physical 
constraints: i) the amount of flatness of the rotational curve of our Galaxy, considering 
conservatively 0.8 • v < v}.^ < 1.2 • v , where v^t is the value of rotational curve at 
distance of 100 kpc from the galactic center; ii) the maximal non dark halo components 
in the Galaxy, considering conservatively 1 • 1O 1O M < M vis < 6 • 1O 1O M [111, 112]. 
In particular, the allowed intervals for p are reported in Table 7. For the models of 
class C and D, as discussed before, only the case of maximal halo (which correspond 
to M vis = 0) has been considered in the analyses discussed in the following. 



7.1.4 The form factors 

In order to take into account the finite dimension of the nucleus in the scattering 
processes, it is necessary to introduce the nuclear form factor, F, that generally also 
depends on the nature of the interaction. 

In the case of SI interactions, we can factorize the total cross section, pointing out 
the contribution given by the form factor: 

cr(q) = asiFljiq), (48) 

here asi is the total cross section when the transfer momentum q is 0. When neglecting 
possible neutron and proton differences, the nuclear form factor can be reasonably 
described by the Fourier transform of charge density p(r) in the nucleus: 



1 f 1 Air f°° 

Fsi(q) = jj P (r)e^ F d 3 r = -— J r S m(qr)p(r)dr. (49) 

In the theoretical calculations and data analysis DAMA/Nal has adopted, for the 
SI form factor, the expression suggested by Helm in [115]: 



c, , x 3ji (qr a ) 

Fsi(q) = exp 

qr 



1 2 2 

"2 S q 



(50) 



where ro = yr^— 5s 2 , r n is the effective nuclear radius, s ~ 1 fm is a parameter 
that allows to take into account the thickness of the nuclear surface and ji(qro) is the 
spherical Bessel function of index 1 . We remind that this expression of SI form factor 
is derived assuming a Fermi distribution for the nuclear charge. 

In Fig. 23 on the right, the effect of a relatively small variation (20%) of the nuclear 
radius, r n , and the nuclear surface thickness parameter, s, in the Helm SI form factor 
is shown; as it can be seen, even a relative small variation of these parameters can 
produce sizeable change in the behaviour of the form factor and, therefore, in the 
expected SI signal rate and in the final result. 

In addition, also other expressions have been considered in literature for the SI 
nuclear form factors of the various nuclei. Just as an example, in Fig. 23 on the 
left, SI form factors discussed e.g. in ref. [116] are depicted for the Iodine nucleus; 
the continuous line represents the nuclear form factor obtained when using expression 
(50) at the fixed assumed values for the related parameters. 
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E R (keVee) E R (keVee) 

Figure 23: Left panel: some SI nuclear form factors considered so far in the literature 
(e.g. in ref. [116]) calculated for the Iodine nucleus. Right panel: example of the 
variation of the SI nuclear form factor, calculated according to ref. [115], when some 
uncertainties on its parameters are considered; as it can be seen, even a relative small 
variation can produce sizeable change in the behaviour of the nuclear form factor 
and, therefore, in the expected SI signal rate. Note that here the notation keVee is 
cxplicitely mentioned in the panels to indicate keV electron equivalent in order to 
remind that the recoil energy has been quenched there by using the quenching factor 
value measured by DAMA in ref. [46] . 



As mentioned, in the analyses given in the following the most cautious Helm SI 
form factor has been adopted taking into account some uncertainties on the nuclear 
radius and on the nuclear surface thickness parameters. As it can be seen, this form 
factor is the less favourable one for Iodine and requires larger SI cross sections for a 
given signal rate. For example, in case the other form factor profiles considered in 
literature would be used, the allowed regions given in the following sections would 
extend to lower cross sections. 

In the case of SD interactions an analytical universal expression for the form factor 
does not exist. In fact, in this case, the internal degrees of the WIMP particle model 
(e.g. supersymmetry in the case of neutralino) cannot be completely decoupled from 
the nuclear ones. Therefore, in order to take into account the property of the nucleus 
interested in the interaction, we have to refer not only to a particular nuclear model 
but also to a particular particle physics model. As an example, if we consider the case 
of the neutralino in the MSSM model, the differential SD cross section can be written 
from eq. (8): 

{^-xy^ mNk2j{j+i)FUEA ,5i) 

where Fsd(Er) is the SD nuclear form factor defined (see e.g. ref. [117]) as 
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f sd( e r) = S(q)/S(0) with ttS(O) = (2 J + 1)A 2 J {J + 1) and: 

S(q) = a 2 S oa (q) + C%S u (q) + a aai S i(q), (52) 

with ao = a p + a n and a\ = a p — a n , where a p and a n are the effective ncutralino- 
nucleon SD coupling strengths defined in §7.1.1. The functions Sij{q) generally depend 
on the considered nuclear model. For nuclei of interest for Dark Matter detection, (as 
the 127 I and 23 Na), the Sij(q) have been parameterized in ref. [117] as function of the 
variable y = (qb/2) 2 (where & is a parameter of the theoretical model, b ~ A 1 / 6 fm): 

8 

S ij (q)=e- 2 y^Ciy l , (53) 

1=0 

where the coefficients C; generally depend on the adopted nucleon-nucleon potential. 
For the 23 Na nucleus only the first three orders of the previous expression have been 
considered in ref. [117], while for the 127 I nucleus two different nuclear potentials (the 
Nijmegen II and the Bonn-A ones) have been used in the evaluation of ref. [117]. 



9=2.435 rad, a n /a p =-0.85 




5 10 15 20 25 30 5 10 15 20 25 30 

E R (keVee) E R (keVee) 

Figure 24: Left panel: some SD nuclear form factors considered so far in the literature 
(e.g. in ref. [116]) calculated for the Iodine nucleus. Right panel: example of variation 
of the Iodine SD nuclear form factor, calculated according to ref. [117] for two different 
choices of the nucleon-nucleon potential and including also a possible small variation 
(20%) of its b parameter. As it can be seen, even a relative small variation either of 
this parameter or of the nuclear potential can produce sizeable change in the behaviour 
of the nuclear form factor and, therefore, in the expected SD signal rate. Note that 
here the notation keVee is explicitely mentioned in the panels to indicate keV electron 
equivalent in order to remind that the recoil energy has been quenched there by using 
the quenching factor value measured by DAMA in ref. [46] . 

It is worth to note that the SD form factor depends on the nature of the interacting 
WIMP and on the nuclear potential. Therefore, the SD form factor is an important 
source of uncertainties in the calculation of the expected rate. As an example, in 
Fig. 24 on the right the effect of the different choice of the nucleon-nucleon potential 
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and of a relatively small variation (20%) of the b parameter in the Iodine SD nuclear 
form factor, calculated according to ref. [117], is shown (we report here, for simplicity, 
only the case with 9 = 2.435 or a n /a p = —0.85, that is pure Zq coupling). Similar 
uncertainties are present for every nucleus. 

There are other expressions considered in literature for the SD nuclear form factor; 
they can be very different. Just as an example Fig. 24 on the left shows some SD form 
factors discussed e.g. in ref. [116] calculated for the Iodine nucleus. 

In the analyses, presented in the following , the SD form factors in the neutralino 
case for Sodium and Iodine nuclei calculated by [117] using - for the Iodine case - 
the nuclear potential by Nijmegen II have been adopted. Analogously, as for the SI 
case, here uncertainties on the b parameter have been included in the evaluations of 
the results. 

For the sake of completeness, let us remind that typically only purely SI or purely 
SD Zq coupling WIMP interactions are considered in the evaluations of the results in 
this field among all the wide available possibilities and that the uncertainty on the 
form factors also largely affects comparisons among results obtained by using different 
target nuclei. 

7.1.5 The spin factors 

Further significant uncertainties in the evaluation of the SD interaction rate arise also 
from the adopted spin factor for the single target-nucleus [117, 118]. As an example, 
the spin factors of some target-nuclei calculated in different models are reported in 
Table 8. Moreover, also for a fixed nuclear model, differences arise from the use - in 
the calculations - of different nuclear potentials as it is e.g. the case of the 129 Xe and 
of the 127 I nuclei, for which similar calculations are already available. 

As it can be noted in Table 8, since the spin factor is a multiplicative factor in the 
expected SD signal rate, its value can also drastically affect the expectations in direct 
search experiments and, therefore, also the inferred exclusion plots or allowed regions 
can largely vary as well as the results of any comparison. 

Moreover, for a complete analysis of a SD component it is worth to remind that 
9 (whose tangent is the ratio between the SD WIMP-neutron and SD WIMP-proton 
effective strengths; see §7.1.1), can continuously assume values in the range to tt. 
For example, in Table 9 spin factors calculated on the basis of ref. [117] are given for 
some 9 values considering few target nuclei and two different nuclear potentials. 

In conclusion, not only large differences in the rate can be expected when using 
target nuclei sensitive to the SD component of the interaction (such as e.g. 23 Na and 
127 1) with respect to those largely insensitive to such a coupling (such as e.g. nat Ge 
and nat Si), but also when using different target nuclei although all - in principle - 
sensitive to such a coupling (compare e.g. the Xenon and Tellurium case with the 
Sodium and Iodine case in Table 9). 

Moreover, other nuclear models and calculations beyond those reported here arc 
possible, introducing large uncertainties in the right estimate of the used spin factor 
for each given target-nucleus. 
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Table 8: Some spin factors estimates assuming simple different models. The values 
given in this Table are A 2 J(J + I) /a 2 ., where a x is either a p or a n depending on the 
unpaired nuclcon. 



Target-Nucleus 


single particle 


odd group 


Comment 


*»Si 


0.750 


0.063 


Neutron is 


73 Ge 


0.306 


0.065 


the unpaired 


129 Xc 


0.750 


0.124 


nucleon 


131 Xc 


0.150 


0.055 






0.750 


0.750 




19 F 


0.750 


0.647 




23 Na 


0.350 


0.041 


Proton is 


27 Al 


0.350 


0.087 


the unpaired 


69 Ga 


0.417 


0.021 


nucleon 


71 Ga 


0.417 


0.089 




75 As 


0.417 


0.000 




127J 


0.250 


0.023 





Table 9: Spin factors calculated on the basis of rcf. [117] for some of the possible 9 
values considering some target nuclei and two different nuclear potentials. The values 
given in this Table are A 2 J(J + l)/a 2 , where a 2 has been defined in §7.1.1. 



Target-Nucleus / 
nuclear potential 


6>=0 


6=7r/4 


6>=7r/2 


61=2.435 
(pure Z 
coupling) 


23 Na 


0.102 


0.060 


0.001 


0.051 


127 I/Bonn A 


0.134 


0.103 


0.008 


0.049 


127 I/Nijmegen II 


0.175 


0.122 


0.006 


0.073 


129 Xe/Bonn A 


0.002 


0.225 


0.387 


0.135 


129 Xe/Nijmegen II 


0.001 


0.145 


0.270 


0.103 


131 Xe/Bonn A 


0.000 


0.046 


0.086 


0.033 


131 Xe/Nijmcgcn II 


0.000 


0.044 


0.078 


0.029 


125 Te/Bonn A 


0.000 


0.124 


0.247 


0.103 


125 Tc/Nijmcgcn II 


0.000 


0.156 


0.313 


0.132 



7.1.6 The quenching factors 

The proper knowledge of other quantities is also necessary for a WIMP direct search 
such as e.g. the recoil/electron response ratio for the given nucleus in the given detector 
(named quenching factor, q). The recoil/electron response ratio can be measured 
with a neutron source or at a neutron generator (see as an example ref. [72]). In 
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the latter case a set-up similar to that reported in Fig. 25 can be used on a quasi- 
monochromatic neutron beam tagging the scattered neutrons. Let us comment that 




Figure 25: Example of a set-up for the measurement of the recoil/electron response 
ratio of a detector. This example is taken from ref. [72]. 

while this configuration works out well for the measurements of recoil/electron response 
ratio, it is unsuitable to determine the absolute detection efficiency for recoils because 
of the uncertain knowledge of the elastic cross sections, of the background due to 
inelastic scatterings, of the electromagnetic and hadronic background from neutron 
interaction in the environment, of the cuts for noise rejection in a very high rate 
environment, of the duty cycle (which is very small), etc. 

Of course, significant differences are often present in literature for the measured 
value of this recoil / electron response ratio for the same nucleus in similar detectors as 
it can be clearly deduced from Table 10, where the quenching factors measured for the 
detectors most commonly used in this field are reported. As it can be seen, significant 
differences in the measured values are present also for the same nucleus in the same 
kind of detectors. This is generally due to different peculiarities of the detectors 
themselves, besides possible additional experimental uncertainties. For example, in 
doped scintillator it can depend on the dopant concentration, in liquid Xenon on the 
residual trace contaminants due to specific experimental features (such as the initial 
purity of the used Xenon gas, the inner surface treatment, the level of vacuum reached 
before filling, the used purification line components and the degassing/release features 
of all the materials of the inner vessel) and in Ge or Si on impurities, etc. Moreover, 
some dependence of the recoil / electron response ratio on the energy has been quoted 
in several cases (see Table 10). 

As far as regards the bolometers, no direct measurement of the recoil energy of 
the target-nucleus has been reported up to now by any of the groups involved in this 
activity, although several bolometers have been irradiated with neutrons along the past 
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Tabic 10: Quenching factor, q, values measured by using neutron sources or neutron 
beams for some detectors and nuclei. When a significant q dependence on the energy 
has been explicitely given, only the maximum and the minimum central values are 
reported for the given energy interval without quoting here - for simplicity - the 
associated errors. 



Nucleus / Detector 


Recoil Energy (keV) 


q 


Reference 


Nal(Tl) 


(6.5-97) 


(0.30 ± 0.01) for Na 


[46] 




(22-330) 


(0.09 ± 0.01) for I 


[46] 




(20-80) 


(0.25 ± 0.03) for Na 


[119] 




(40-100) 


(0.08 ± 0.02) for I 


[119] 




(4-252) 


(0.275 ± 0.018) for Na 


[120] 




(10-71) 


(0.086 ± 0.007) for I 


[120] 




(5-100) 


(0.4 ± 0.2) for Na 


[121] 




(40-300) 


(0.05 ± 0.02) for I 


[121] 


CaF 2 (Eu) 


(30-100) 


(0.06-0.11) for Ca 


[120] 




(10-100) 


(0.08-0.17) for F 


[120] 




(90-130) 


(0.049 ± 0.005) for Ca 


[45] 




(75-270) 


(0.069 ± 0.005) for F 


[45] 




(53-192) 


(0.11-0.20) for F 


[122] 




(25-91) 


(0.09-0.23) for Ca 


[122] 


CsI(Tl) 


(25-150) 


(0.15-0.07) 


[123] 




(10-65) 


(0.17-0.12) 


[124] 




(10-65) 


(0.22-0.12) 


[125] 


CsI(Na) 


(10-40) 


(0.10-0.07) 


[125] 


Ge 


(3-18) 


(0.29-0.23) 


[126] 




(21-50) 


(0.14-0.24) 


[127] 




(10-80) 


(0.18-0.34) 


[128] 




(20-70) 


(0.24-0.33) 


[129] 


Si 


(5-22) 


(0.23-0.42) 


[130] 




22 


(0.32 ± 0.10) 


[131] 


Liquid Xe 


(30-70) 


( 0.46 ± 0.10) 


[72] 




(40-70) 


(0.18 ± 0.03) 


[132] 




(40-70) 


(0.22 ± 0.01) 


[133] 


Bolometers 




none available 








at time of writing this paper 








(assumed 1) 





decade. For the sake of completeness, we remind that a measurement of the response 
of a Tc02 bolometer to surface 224 Ra recoiling nuclei has been reported in rcf. [134]; 
this measurement, although its importance, does not represent a determination of the 
quenching factor of the target- (either Tc or O) nuclei of the Tc02 bolometer. In fact, 
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the recoiling nuclei arc not the Tc and O ones and the "external" recoils are generated 
on the detector surface where the sensors are located and, thus, do not involve the 
response of the whole bulk of the target-detector. Anyhow, these values cannot of 
course be extended to whatever kind of bolometer. 

7.1.7 Some miscellaneous 

Besides the uncertainties already discussed, there exist a large number of experimental 
details that have to be considered in the calculations of the signal expectations. For 
example, it must be suitably studied, checked, monitored and discussed the role played 
by the stability of the energy scale, by its right evaluation (see also §7.1.6), by a reliable 
identification of the energy threshold and of the residual noise above it (as it can be 
effectively done in Nal(Tl) detectors with adequate number of photoelectrons/keV, 
see e.g. [41]), by external veto (especially when "high" rate anticoincidences are used) 
and by the rejection procedures used in some experiments to filter the data, etc.. In 
the case of DAMA/Nal details have been given along the last about ten years and 
published (see e.g. [41, 61, 63, 64]). 

Finally, let us remind that - when results obtained by using different target nuclei 
are considered - also the effect of the uncertainty on the scaling laws of the nuclear 
cross sections to the reference ones (e.g., generally, on nucleon) must be taken into 
account. In fact, it is common practice to use the WIMP-nucleon cross sections a si 
and usd and the scaling laws reported in §7.1.1; but - in principle - other scaling laws 
cannot be excluded at the present knowledge of the real nature of a WIMP candidate. 
This can be an additional uncertainty in quests for a candidate and in comparisons 
among experiments using different target nuclei. 

7.1.8 Priors 

It is common practice in extracting physical information from the data to account for 
related priors. 

In particular, in the quest for the candidate particle two main priors have been 
considered in the DAMA/Nal first quests for a candidate [60, 61, 62, 63, 65, 66, 11]. 

The first prior, which has been properly included in rcf. [63, 65, 66, 11] and 
also considered in the following, accounts for the upper limits measured on the recoil 
fractions in the data of the DAMA/NaI-0 running period [46, 64], which was carefully 
and especially devoted to such an investigation. 

The second prior regards the mass limit for supersymmetric candidates, achieved 
- within some assumed model frameworks - by experiments at accelerators. In par- 
ticular, because of this prior, WIMP masses above 30 GeV (25 GeV in ref. [60]) have 
been investigated in refs. [61, 63, 65, 66, 11] for few (of the many possible) model 
frameworks. Specifically, it accounted for the lower bound on the neutralino mass as 
derived from the LEP data in the adopted supersymmetric schemes based on GUT 
assumptions [78]. However, other model assumptions are possible and would imply 
significant variations of some accelerators bounds. As an example, we mention the 
recent ref. [135] where the assumption on the gaugino-mass unification at GUT scale 
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has been released 13 . The development of these schemes is very interesting since - as 
well known - DAMA/Nal is intrinsically sensitive both to low and high WIMP mass 
having both a light (the 23 Na) and a heavy (the 127 I) target-nucleus. 

However, still following the present model dependent results quoted by LEP in the 
supersymmetric schemes based on GUT assumptions the considered lower bound is at 
present 37 GeV [138] 14 . It worth to note that this mass limit selects the WIMP-Iodine 
elastic scattering as dominant because of the adopted scaling laws and of kinematical 
arguments. 

7.2 Results on the quest for a candidate in some of the possible 
model frameworks 

Just as a corollary of the model independent result given in §6.1, in the following 
some of the many possible model dependent quests for a WIMP candidate is carried 
out using the data collected during all the seven annual cycles and considering all 
the halo models summarized in §7.1.3 for three of the possible values of the local 
velocity v?,: 170 km/s, 220 km/s and 270 km/s. The used halo density follows the 
prescriptions of §7. 1.3. VI. The escape velocity has been maintained at the fixed value: 
650 km/s; of course, it is worth to note that the present existing uncertainties affecting 
the knowledge of the escape velocity will significantly extend allowed regions e.g. in 
the cases of preferred inelastic WIMPs and of light mass WIMP candidates; its effect 
would be instead marginal at large WIMP masses (see e.g. the case for exclusion plots 
given in Fig. 3). 

In particular, possible scenarios have been exploited for the halo models described 
in §7.1.3 in some discrete cases either considering the mean values of the parameters 
of the used nuclear form factors and of the measured quenching factors (case A) or 
adopting the same procedure as in refs. [65, 66] 15 (case B) or in one of the possible 
more extreme cases where the Iodine nucleus parameters are fixed at the values of case 
B, while for the Sodium nucleus one considers: i) 23 Na quenching factor at the lowest 
value measured in literature (see Table 10); ii) the nuclear radius, r n , and the nuclear 
surface thickness parameter, s, in the SI Form Factor [115] from their central values 
up to +20%; iii) the b parameter in the considered SD form factor from the given value 
[117] up to +20% (case C). 

In the following sections, for simplicity, the results of these corollary quests for a 
candidate particle is presented in terms of allowed regions obtained as superposition 
of the configurations corresponding to likelihood function values distant more than 4cr 
from the null hypothesis (absence of modulation) in each of the several (but still a 

13 In this case also neutralino masses down to ~ 6 GeV are possible, this lower bound being de- 
termined by current upper limit on relic abundance for cold dark matter (somewhat higher values of 
15-18 GeV are obtained for the neutralino mass lower bound, if Higgs masses are assumed large ( 1 
TeV) [136, 137]. 

14 Highcr limits are available for other more constrained models. These latter ones as well as 
possible future increase of the present 37 GeV lower bound would further select - for these scenarios - 
the possible models in the quest for the candidate from the DAMA/Nal data favouring, in particular, 
halo models with small local velocity and/or co-rotation. 

15 that is, by varying either: i) the mean values of the measured 23 Na and 127 I quenching factors 
[46] up to +2 times the errors; ii) the nuclear radius, r„, and the nuclear surface thickness parameter, 
s, in the SI Form Factor [115] from their central values down to -20%; iii) the b parameter in the 
considered SD form factor from the given value [117] down to -20%. 
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limited number) of the possible model frameworks considered here. Priors have been 
discussed in §7.1.8. 

Obviously, larger sensitivities than those reported in the following would be reached 
when including the effect of other existing uncertainties on assumptions and related 
parameters, as it can be also inferred from the previous sections. 

7.2.1 WIMPs with mixed SI&SD interaction in some of the possible model 
frameworks 

The most general scenario of WIMP nucleus elastic interaction, to which the 
DAMA/Nal target nuclei are fully sensitive, is the one where both the SI and the 
SD components of the cross section (see §7.1.1) are present. Thus, as first we intro- 
duce here the case for a candidate with both SI and SD couplings to ordinary matter 
similarly as in rcf. [65]. 

As already described in §7.1.1, in this most general scenario the space of the free 
parameters is a 4-dimensional volume defined by mw, £&si, £,csd and 9 (which varies 
from to 7r). Thus, the general solution would be a four dimensional allowed volume 
for each considered model framework. Since the graphic representation of this allowed 
volume is quite difficult, we show in Fig. 26 the obtained regions in the plane t;crsi vs 
£&sd for some of the possible 9 and mw values in the model frameworks considered 
here. In particular, we report just four couplings, which correspond to the following 
values of the mixing angle 9: i) 9 = (a n =0 and a p ^ or \a p \ >> \a n \) corresponding 
to a particle with null SD coupling to neutron; ii) B = ir/4 (a p = a n ) corresponding to 
a particle with the same SD coupling to neutron and proton; iii) 9 = tt/2 (a n ^ and 
a p = or \a n \ » | a p |) corresponding to a particle with null SD couplings to proton; 
iv) 9 — 2.435 rad (^ = -0.85) corresponding to a particle with SD coupling through 
Z exchange. The case a p — —a n is nearly similar to the case iv). 

To offer an example of how the allowed regions have been built, Fig. 27 shows 
cxplicitely the superposition of the slices obtained for each one of the model frameworks 
considered here in the particular case of mw — 90 GeV and 9 = 2.435 (pure Zq 
coupling). 

From the given figures it is clear that at present either a purely SI or a purely SD 
or a mixed SI&SD configurations are supported by the experimental data of the seven 
annual cycles. 

Some other comments related to the effect of a SD component different from zero 
will be also addressed in the following. 

7.2.2 WIMPs with dominant SI interaction in some of the possible model 
frameworks 

Generally, mainly the case of purely spin-independent coupled WIMP is considered 
in literature. In fact, often the spin-independent interaction with ordinary matter is 
assumed to be dominant since e.g. most of the used target-nuclei are practically not 
sensitive to SD interactions (as on the contrary 23 Na and 127 I are) and the theoretical 
calculations are even much more complex and uncertain. 
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Figure 26: A case of a WIMP with mixed SI&SD interaction in the model frameworks 
given in the text. Colored areas: example of slices (of the allowed volume) in the 
plane £,asi vs £ctsd for some of the possible mw and 9 values. See §7.2. Inclusion 
of other existing uncertainties on parameters and models (as previously discussed to 
some extent in this paper) would further extend the regions; for example, the use of 
more favourable form factors than those we considered here (see §7.1.4) alone would 
move them towards lower cross sections. 



Thus, following an analogous procedure as for the previous case, we have exploited 
for the same model frameworks the purely SI scenario alone. In this case the free 
parameters are two: mw and t;o~si- 

In Fig. 28 the region allowed in the plane mw and £,o~si for the considered model 
frameworks is reported. The vertical dotted line represents the model dependent prior 
discussed in §7.1.8, that is the present lower bound on supersymmetric candidate as 
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^g sd (pb) 

Figure 27: This figure explicitely shows all the slices (of the allowed volume) in the 
plane £05/ vs £<jsd obtained for mw — 90 GeV and 9 = 2.435 (pure Z coupling) 
for each one of the considered model frameworks (see §7.2). The region included in 
between the two extreme lines is that shown for the same mw and 9 values in Fig. 26 
(see also the related caption). 

derived from the LEP data in supersymmetric scheme with gaugino-mass unification 
at GUT (see §7.1.8). The configurations below the vertical line can be of interest 
for neutralino when other schemes are considered (see §7.1.8) and for generic WIMP 
candidate. As shown in Fig. 28, also WIMP masses above 200 GeV are allowed, 
in particular, for every set of parameters' values when considering low local velocity 
and: i) the Evans' logarithmic CI and C2 co-rotating halo models; ii) the triaxial D2 
and D4 non-rotating halo models; iii) the Evans power-law B3 model, but only with 
parameters as in set A). 

Of course, best fit values of cross section and WIMP mass span over a large range 
depending on the model framework. Just as an example, in the triaxial D2 halo model 
with maximal p , v = 170 km/s and parameters as in the case C, the best fit values 
are m w = (74±^) GeV and £a SI = (2.6 ± 0.4) • 10~ 6 pb. 

Effect of a SD component different from zero on allowed SI regions 

Let us now point out, in addition, that configurations with £05/ even much lower 
than those shown in Fig. 28 would be accessible also if an even small SD contribution 
would be present in the interaction as described in §7.2.1. This possibility is clearly 
pointed out in Fig. 29 where an example of regions in the plane (mw, £,&si) corre- 
sponding to different SD contributions are reported for the case 9 = 0. In this example 
the Evans' logarithmic axisymmetric C2 halo model with vq = 170 km/s, po equal to 
the maximum value for this model (see Table 7) and the set of parameters A have 
been considered. The values of £<jsd range there from to 0.08 pb. As it can be seen, 
increasing the SD contribution the regions allowed in the (mw, C^S/) plane involve SI 
cross sections much lower than 1 x 10~ 6 pb. It can be noted that for cjsd > 0.08 pb the 
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200 400 

m w (GeV) 

Figure 28: Case of a WIMP with dominant SI interaction for the model frameworks 
given in the text. Region allowed in the plane (mw, C^s/)- See §7.2; the vertical 
clotted line represents the model dependent prior discussed in §7.1.8. The area at 
WIMP masses above 200 GeV is allowed for low local velocity - wo=170km/s - and all 
considered sets of parameters by the Evans' logarithmic CI co-rotating halo model, 
by the Evans' logarithmic C2 co-rotating halo model, by the triaxial D2 and DA non- 
rotating halo models and also by the Evans power-law B3 model with parameters of 
the set A). The inclusion of other existing uncertainties on parameters and models (as 
previously discussed to some extent in this paper) would further extend the region; for 
example, the use of more favourable SI form factor for Iodine (see §7.1.4) alone would 
move it towards lower cross sections. 

annual modulation effect observed is also compatible - for mw — 40 — 75 GeV - with 
a WIMP candidate with no SI interaction at all in this particular model framework. 

These arguments clearly show that also a relatively small SD contribution can 
drastically change the allowed region in the (mw, £,o~Si) plane; therefore, e.g. there is 
not meaning in the bare comparison between regions allowed in experiments that are 
also sensitive to SD coupling and exclusion plots achieved by experiments that are not. 
The same is when comparing regions allowed by experiments whose target-nuclei have 
unpaired proton with exclusion plots quoted by experiments using target-nuclei with 
unpaired neutron when the SD component of the WIMP interaction would correspond 
either to 9 ~ or 9 ~ n. 

7.2.3 WIMPs with dominant SD interaction in some of the possible model 
frameworks 

Let us now focus on the case of a candidate with purely spin-dependent coupling to 
which DAMA/Nal is - as mentioned - fully sensitive. 

When the SD component is different from zero, a very large number of possible 
configurations is available (see §7.1.1). In fact, in this scenario the space of free pa- 
rameters is a 3-dimensional volume defined by mw, £,&sd and 9 (which can vary from 
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Figure 29: Example of the effect induced by the inclusion of a SD component different 
from zero on the allowed regions given in the plane £<isi vs mw In this example the 
Evans' logarithmic axisymmctric C2 halo model with vq = 170 km/s, po equal to the 
maximum value for this model (see Table 7) and the set of parameters A have been 
considered. The different regions refer to different SD contributions for the particular 
case of 6 = 0: a S D = pb (a), 0.02 pb (b), 0.04 pb (c), 0.05 pb (d), 0.06 pb (e), 
0.08 pb (f). See §7.2; the vertical dotted line represents the model dependent prior 
discussed in §7.1.8. 



to 7r). Here, for simplicity as already done in §7.2.1, we show the results obtained 
only for 4 particular couplings, which correspond to the following values of the mixing 
angle 9: i) 9 = (a n —0 and a p ^ or \a p \ >> \a n \; ii) 9 = 7r/4 (a p = a n ); iii) 9 = 
tt/2 (a n ^ and a p = or \a n \ » \a p \; iv) 9 = 2.435 rad = -0.85). 

Fig. 30 shows the regions allowed in the plane (mw, &sd) for the same model 
frameworks quoted above; other configurations are possible varying the 9 value. The 
area at WIMP masses above 200 GeV is allowed for low local velocity - t>o = 170km/s 
- and all considered sets of parameters by the Evans' logarithmic C2 co-rotating halo 
model. 

Moreover, the accounting for the uncertainties e.g. on the spin factors as well as 
different possible formulations of the SD form factors would extend the allowed regions, 
e.g. towards lower £<jsd values. 

Finally, (,<7sd lower than those corresponding to the regions shown in Fig. 30 are 
possible also e.g. in case of an even small SI contribution, as shown in Fig. 31. 



7.2.4 WIMPs with preferred inelastic interaction in some of the possible 
model frameworks 

An analysis considering the same model frameworks has been carried out for the case 
of WIMPs with preferred inelastic interaction (see §7.1.2). 

In this inelastic Dark Matter scenario an allowed volume in the space (£a p ,mw,8) 
is obtained. For simplicity, Fig. 32 shows slices of such an allowed volume at some 
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Figure 30: Case of a WIMP with dominant SD interaction in the model frameworks 
given in the text. Regions allowed in the plane (mw, £,o~sd)- See §7.2; the vertical 
dotted line represents the model dependent prior discussed in §7.1.8. The panels refer 
to only few particular cases for 9 (which can instead vary between and 7r). The 
area at WIMP masses above 200 GeV is allowed for low local velocity - t>o=170km/s 
- and all considered sets of parameters by the Evans' logarithmic C2 co-rotating 
halo model. Inclusion of other existing uncertainties on parameters and models (as 
previously discussed to some extent in this paper) would further extend the regions; 
for example, the use of more favourable SD form factors (see §7.1.4) alone would move 
them towards lower cross sections. 



given WIMP masses. 

There the superpositions of the allowed regions obtained, when varying the model 
framework within the considered set, are shown for each mw- As a consequence, 
the cross section value at given S can span over several orders of magnitude. The 
upper border of each region is reached when vthr approximates the maximum WIMP 
velocity in the Earth frame for each considered model framework. It can also be 
noted that when mw 3> m^, the expected differential energy spectrum is trivially 
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Figure 31: Example of the effect induced by the inclusion of a SI component different 
from zero on the allowed regions in the plane £osd vs mw In this example the 
Evans' logarithmic axisymmetric C2 halo model with vq = 170 km/s, pa equal to the 
maximum value for this model (see Table 7) and the set of parameters A for 9 = 
have been considered. The different regions refer to different SI contributions with: 
a SI = pb (a), 2 x 10~ 7 pb (b), 4 x 10~ 7 pb (c), 6 x 10~ 7 pb (d), 8 x 10~ 7 pb (e), 
10~ 6 pb (f). Sec §7.2; the vertical dotted line represents the model dependent prior 
discussed in 57. 1.8. 



dependent on mw and, in particular, it is proportional to the ratio between £a p and 
myy] therefore for very high mass the allowed region can be obtained straightforward. 
We remind that in these calculations v esc has been assumed at fixed value, while its 
present uncertainties can play a significant role in the scenario of WIMP with preferred 
inelastic scattering as mentioned in §7.1.2. 

Note that each set of values (within those allowed by the associated uncertainties) 
for the previously mentioned parameters gives rise to a different expectation, thus to 
a different best fit values. As an example we mention the best fit values for mw = 70 
GeV in the NFW B5 halo model with vo = 170 km/s, maximal po in this model and 
parameters as in case B): (5 = 86tg) keV and £a p = (1.2 ± 0.2) x 10~ 5 pb. 

7.2.5 Conclusion on the quest for a candidate in some of the possible 
model frameworks 

In this section the possible nature of a candidate, which could account for the observed 
model independent evidence, has been investigated by exploring - as already done on 
the partial statistics [60, 61, 62, 63, 64, 65, 66, 11] - various kinds of possible couplings 
and some (of the many) possible model frameworks. 

We stress that, although several scenarios have been investigated, the analyses are 
not exhaustive at all of the existing possibilities because of the poor present knowledge 
on many astrophysical, nuclear and particle physics assumptions and related parame- 
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Figure 32: Case of a WIMP with preferred inelastic interaction in the model frame- 
works given in the text. Examples of slices (colored areas) of the allowed volumes (^cr p , 
S, mw) for some m-w values for the considered model frameworks for WIMP with pre- 
ferred inelastic interaction. See §7.2. In these calculations v esc has been assumed at 
fixed value, while its present uncertainties would play a significant role in the obtained 
results. Inclusion of other existing uncertainties on parameters and models (as pre- 
viously discussed to some extent in this paper) would further extend the regions; for 
example, already the use of a more favourable SI form factor for Iodine (see §7.1.4) 
alone would move them towards lower cross sections. 



ters as well as of the existing uncertainties in the determination of some experimental 
parameters which are necessary in the calculations. For example, other parameters 
values can be considered for the investigated halo models as well as other different 
halo models too, other form factors and related parameters, other spin factors etc.. 
We remind that analogous uncertainties are present in every model dependent result 
(such as e.g. exclusion plots and WIMP parameters from indirect searches); thus, 
intrinsically, bare comparisons have always only a very relative meaning. 

The discussion, carried out in this section, has also allowed to introduce the 
main general arguments related to the model dependent calculations in WIMP di- 
rect searches. 
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8 Conclusion 



In this paper general aspects of the Dark Matter direct search have been reviewed in 
the light of the activity and results achieved by the DAMA/Nal experiment at the 
Gran Sasso National Laboratory of I.N.F.N.. DAMA/Nal has been a pioneer experi- 
ment running at LNGS for several years and investigating as first the WIMP annual 
modulation signature with suitable sensitivity and control of the running parameters. 
During seven independent experiments of one year each one, it has pointed out the 
presence of a modulation satisfying the many peculiarities of a WIMP induced effect, 
reaching a significant evidence. As a corollary result, it has also pointed out the com- 
plexity of the quest for a WIMP candidate because of the present poor knowledge on 
the many astrophysical, nuclear and particle physics aspects. 

As regards other experiments - to have a realistic comparison - experiments in- 
vestigating with the same sensitivity and control of the running condition the annual 
modulation signature are necessary. Of course, the target nuclei also play a crucial role, 
since they can offer significantly different sensitivities depending e.g. on the nature of 
the WIMP particle and on their nuclear properties. 

The growing in the field of serious and independent efforts searching for WIMP 
model independent signatures will certainly contribute to increase the knowledge in 
the field as well as efforts to more deeply investigate models and parameters. 

Some of the most competitive activities for the near future, exposing a significantly 
large target-mass, are starting at the Gran Sasso National Laboratory: CUORICINO, 
GENIUS-TF (which will also be devoted to the investigation of double beta decay 
processes) and our new experiment DAMA/LIBRA. In fact, on our behalf, after the 
completion of the data taking of the ~ 100 kg Nal(Tl) set-up (on July 2002), as a 
result of our continuous efforts toward the creation of ultimate radiopure set-ups, the 
new DAMA/LIBRA has been installed (see ref. [139]). The LIBRA set-up is made by 
25 Nal(Tl) detectors, 9.70 kg each one. The new detectors have been realised thanks 
to a second generation R&D with Crismatec/SaintGobain company, by exploiting in 
particular new radiopurification techniques of the Nal and Til selected powders. In the 
framework of this R&D new materials have been selected, prototypes have been built 
and devoted protocols have been fixed and used. The whole installation has largely 
been modified. This new DAMA/LIBRA set-up, having a larger exposed mass and an 
higher overall radiopurity, will offer a significantly increased sensitivity to contribute 
to further efforts in improving the understanding of this field. 
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